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ABSTRACT 

29       30 
The  excitation  function  for  the  Si   (p,  ft)P   reaction  has  been 

experimentally  determined  in  the  proton  energy  range  from  1420  kev 
through  2160  kev.   Unreported  resonances  in  this  reaction  have  been 
identified  at  1975,  2033,  2075,  and  2117  kev  with  estimated  errors 
of  +  4  kev.   Gamma  spectra  at  these  resonances  are  obscured  by 
contaminant  radiation  but  suggest  complex  decay  schemes.   Double 
angular  correlation  measurements  have  been  made  at  the  1470-kev  re- 
sonance indicating  a  resonance  level  spin  of  2+,  or  less  probably, 
2~.   The  level  at  4.50  Mev  is  shown  to  probably  have  a  spin  of  2+, 
or  less  probably,  3+  rather  than  0+  as  has  been  predicted. 
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I.   INTRODUCTION 

As  a  result  of  successive  refinements  and  extensions  of  the  collec- 
tive model  of  nuclear  structure,  it  has  been  possible  to  extend  its  ap- 
plicability into  the  region  of  intermediate  mass  nuclei  (20^A-440)„   By 

considering  the  outer  nucleons  to  be  strongly  coupled  to  the  nuclear  core, 

25    25    28        29 
excited  states  in  Al   ,  Mg   ,  Al   ,  and  Si   have  been  interpreted  as  col- 

1-4 
lective  rotational  excited  states  with  considerable  success.     More  re- 

31       30 
cently  states  in  P   and  Si   have  been  identified  with  collective  vi- 

brational  states  of  a  weakly  coupled  collective  model.  '    The  experimental 
collection  of  information  on  the  states  of  intermediate  mass  nuclei  serves 
to  verify  the  theoretical  models,  to  determine  the  parameters  which  enter 
into  descriptions  of  the  models,  and  to  permit  the  extensions  and  refine- 
ments to  theory  necessary  for  a  more  adequate  description  of  nuclei.   Of 
particular  relevance  to  the  present  investigation  is  the  prediction  based 

on  the  weakly  coupled  collective  model  applied  to  A  ~  30  nuclei  that  the 

30  6 

4.50-Mev  level  of  P   has  a  spin-parity  of  0+.    In  the  present  work  the 

experimental  determination  of  this  value  has  been  the  object  of  double 

angular  correlation  measurements  performed  on  the  gamma  rays  observed  at 

29       30 
the  1470-kev  resonance  in  the  Si   (p,2f)P   reaction. 

29       30 
Determinations  of  resonance  energies  in  the  Si   (p3ft)P   reaction 

have  been  carried  out  by  numerous  investigators  since  1941.  The  pro- 
ton energies  which  have  been  employed  in  these  investigations  extend  over 
the  range  300  to  1950  kev. 

In  addition  to  detecting  resonances  and  measuring  the  energies  at 
which  they  occur,  the  allied  problem  of  analyzing  the  associated  gamma-ray 
spectra  has  been  the  object  of  an  increasingly  large  amount  of  research 
commensurate  with  improvements  in  gamma-ray  spectroscopy  techniques  and 
instrumentation. 

Investigations  of  the  gamma  rays  which  arise  from  deexcitation  of 

30 
the  excited  states  of  P   were  begun  in  1954  by  Endt,  Kluyver,  and  Van 

-i  /*  0 Q  on 

der  Leun.    They  studied  the  Si   (p,tf)P   reaction  at  resonances  oc- 
curring at  proton  energies  of  414  kev  and  326  kev  and  assigned  spin, 


parity,  and  isotopic  spin  to  the  ground  state  and  the  first  excited  state 

17 


30 

of  P   premised  on  gamma- transition  probabilities,  shell  model  consider- 


ations and  the  isotopic  spin  multiplet  formalism.   In  1956  Broude  et  al. 


confirmed  the  assignments  to  the  first  excited  level  by  measuring  the 

angular  distribution  of  the  gamma-rays  at  the  414-kev  resonance.   They 

30 
also  detected  five  new  states  in  P   by  examining  gamma  radiation 

arising  from  the  same  reaction  at  resonances  occurring  at  proton  ener~ 

gies  of  696  kev  and  737  kev. 

1 7  18 

Investigations  by  Lee  and  Mooring  '  and  by  Endt  and  Paris   ,  who 

32      30 
studied  the  S   (d,a)P   reaction^  showed  the  first  excited  state  to  be 

in  reality  a  close  doublet.   The  latter  investigators  were  able  to  accu- 

30 
rately  determine  the  energies  of  30  excited  levels  of  P   up  to  excitation 

energies  of  5.8  Mev  by  magnetic  analysis  of  the  alpha-particle  groups  from 

this  reaction . 

19 
In  1958  Van  der  Leun  and  Endt   reinvestigated  the  four  aforemen- 
tioned resonances,  studying  the  decay  schemes  by  analysis  of  gamma-ray 
spectra,  coincidence  measurements ,  and  angular  distribution  measurements. 
From  this  work  came  decay  schemes,  branching  ratioss  and  assignment  of 

spin-parities  and  isotopic  spins  to  the  four  resonance  levels  and  to  five 

30 
other  excited  levels  of  P   s  including  the  two  states  of  the  above  men- 
tioned doublet.   During  1962  the  gamma-ray  spectra  associated  with  the 

20 
proton  resonances  at  1330  kev   and  the  resonances  at  1375  kev  and  1505 

21 
kev  were  investigated „    Proposed  decay  schemes  were  advanced  for  these 

30 
resonances  and  spin-parity  assignments  made  to  three  levels  of  P   .   One 

29       30 
of  the  more  comprehensive  investigations  of  the  Si   (Ps/OP   reaction 

22 
is  that  of  Baart  et  al.   ,  who  studied  six  resonances  and  proposed  decay 

schemes  for  them.   Firm  spin-parities  were  assigned  to  six  low-  lying 

levels  and  tentative  assignments  were  made  to  four  more  levels.   The 

23 
most  recent  investigation  of  the  reaction  is  that  by  Moore.    As  a  re- 
sult of  this  work  decay  schemes  were  proposed  for  resonances  at  1302, 
1470,  1505,  1640,  16863  1748,,  and  1772  kev,  and  spin-parities  were  de- 
termined for  three  resonance  levels  and  two  intermediate  levels.   The 

30 
present  knowledge  of  the  quantum  states  of  the  lower  P   levels  re- 
sulting from  the  above  investigations  is  summarized  in  Fig.  1.   Two  com- 
pilations which  summarize  important  experimental  results  of  the  above 

24 
mentioned  investigations  are  the  Nuclear  Data  Sheets   and  the  compi- 
lation of  information  on  the  energy  levels  of  nuclei  with  Z  =  11  to 
Z  =  20  by  Endt  and  Van  der  Leun.  5 


Energies  in  MEV,  T  values  are  isotopic  spins, 
spin-parities  in  parentheses  are  less  likely  alternates. 
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THE  LOWER  ENERGY  LEVELS    OF   P 
FIGURE      I. 
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The  purpose  of  the  present  work  has  been  to  investigate  resonances 

29       30 
in  the  Si  (p,1f)¥       reaction  in  the  proton  energy  range  from  1950  kev 

to  2150  kev  which  has  not  previously  been  reported  on,  and  to  investigate 

the  ganuna-ray  spectrum  associated  with  the  known  resonance  occurring  at 

1470  kev.   In  particular,  the  angular  distributions  of  gamma  rays  at  the 

latter  resonance  have  been  measured  in  a  double  angular  correlation  ex- 

30 
periment  to  obtain  information  on  the  4.50-Mev  level  in  P   .   Additionally, 

some  information  has  been  obtained  from  examination  of  spectra  taken  at  the 

19 
higher  resonances  although  contamination  by  F   hampered  this  work.   The 

following  sections  of  this  paper  describe  the  experimental  equipment,  pro- 
cedure, and  data  analysis  employed,  and  the  results  which  have  been  ob= 
tained  from  the  investigation. 


II.      DESCRIPTION  OF   EQUIPMENT 

1.  Experimental  Layout 

Fig.  2  indicates  the  general  arrangement  of  equipment  used  in  this 
experiment.  Fig.  3  is  a  photograph  of  the  target  area  showing  the  cor- 
relation table  tilted  into  a  vertical  plane. 

2.  Accelerator 

The  2-Mev  Van  de  Graaff  horizontal  electrostatic  accelerator  of  the 
U.S.  Naval  Postgraduate  School,  manufactered  by  the  High  Voltage  Engi- 
neering Corporation,  was  used  in  this  work  as  a  source  of  protons.   The 

9  ft 
radio  frequency  ion  source  is  of  the  type  described  by  Moak.    A  proton 

spin-resonance  fluxmeter  for  accurately  measuring  the  analyzing  magnet's 

field  was  installed  furing  the  course  of  the  work.   An  error  signal  from 

the  fluxmeter  was  used  in  a  feed  back  loop  to  the  power  supply  of  the 

beam  analyzing  magnet  to  hold  the  field  strength  automatically  at  a  set 

value.   This  feature  simplified  the  maintenance  of  a  given  resonance 

energy  for  long  bombardment  times  without  "slipping  off", 

A  current  integrator  of  the  vacuum  tube  electrometer  type  was  used 

to  measure  beam  current.   High  stability  plastic  insulated  capacitors 

were  used  to  collect  the  charge.   The  voltage  across  the  capacitor  in 

use  was  applied  to  the  grid  of  a  954  acorn  tube.,  whose  output  fired  a 

Schmidt  trigger  circuit  at  a  preselected  voltage  level.   The  trigger 

output  actuated  other  circuitry  to  discharge  the  capacitor  and  advance 

a  Veeder-Root  cycle  counter  after  each  cycle. 

3.  Target  Assemblies 

Two  types  of  targets  and  target  mounts  were  developed  and  used  in 
this  experiment o   These  are  shown  in  Fig.  4.   For  resonance  and  coinci- 
dence measurements ,  which  were  made  in  the  plane  perpendicular  to  the 
beam,  the  target  backing  consisted  of  a  one-half -in0-diameter  flanged 
copper  hemisphere  held  against  a  flat  copper  boss  on  the  mount  by  a 
threaded  retaining  ring.   This  assembly  was  designed  with  the  beam  di- 
rection as  a  symmetry  axis  so  that  it  can  be  used  for  triple  correlation 
experiments  in  the  perpendicular  plane  with  no  perturbations  from  target 
geometry.   For  the  double  correlations s  measured  in  the  horizontal  plane 
containing  the  beam,,  target  backings  were  0.046  cm  thick  copper  strips 
clamped  in  a  copper  mount.   Both  types  of  backings  and  mounts  were  vacuum 
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plated  with  a  layer  of  gold  thick  enough  to  stop  incident  protons.   For 
both  geometries  the  target  mount  was  fitted  inside  one  end  of  a  one~in. 
diameter  tubular  chamber,,  the  axis  of  which  coincided  with  the  beam  path. 
Both  target  assemblies  were  designed  so  as  to  keep  mass  in  the  target 
vicinity  to  a  minimum.   Competing  with  this  requirement  was  the  need  to 
cool  the  target  in  order  to  prevent  any  appreciable  target  material  eva- 
poration.  Maximum  possible  beam  intensities  were  desired  so  as  to  take 
advantage  of  the  better  counting  statistics  which  accompany  high  yields. 
Still  another  requirement  of  angular  correlation  measurements  is  a  small 
target  spot  to  approximate  a  point  source  of  gamma  rays.   A  beam  colli- 
mator consisting  of  two  five-mil  tantalum  discs  separated  by  0.4  in.  was 
mounted  in  a  brass  holder  38  in.  from  the  target.   A  0.096-in.  diameter 
collimating  aperture  was  followed  by  a  0.098-in.  hole  in  the  second  disc 
for  the  purpose  of  stopping  protons  scattered  from  the  edge  of  the  col- 
limating disc.   A  cylindrical  lead  shield  surrounded  the  collimator  in 
order  to  attenuate  any  gamma  rays  produced  hare. 

As  is  described  in  Chapter  ill;,  target  material  was  evaporated  di- 
rectly onto  the  gold-plated  copper  backings  to  permit  good  heat  conduction, 
Initially,  a  hollow  copper  target  backing  through  which  coolant  could  be 
circulated  was  tried,  but  was  abandoned  for  a  simpler  solid  copper  back- 
ing cooled  by  conduction  through  a  finned  copper  extension  leading  to  an 
external  Dewar  flask.   Liquid  air  in  the  Dewar  proved  to  be  convenient 
as  a  coolant  at  high  beam  powers  and3  below  about  ten  watts,  target  power 
dissipation  by  forced  air  cooling  of  the  target  extension  fins  proved 
adequate. 

After  unsuccessful  trials  of  Lucite  and  glass  target  tubes,  aluminum 
was  finally  adopted  for  this  purpose.   An  auxiliary  target  assembly  with 
a  quartz  window  target  was  used  to  accurately  center  the  proton  beam  spot. 
The  target  chamber  tube,  which  was  insulated  from  ground  at  both  ends., 
was  maintained  at  a  negative  potential  of  300  volts  relative  to  ground  in 
order  to  suppress  secondary  electrons  from  the  target.   In  order  to  pre- 
vent recoil  protons  from  exciting  aluminum  resonances,  the  target  tube 
was  lined  with  one-mil  gold  foil. 

The  hemispherical  copper  target  backings  for  measurements  in  the 
plane  perpendicular  to  the  beam  were  needed  in  some  quantity,  so  were 
"mass  produced"  in  a  small  hand  forge  made  for  this  purpose.   A  surface 
which  needed  only  minor  polishing  and  cleaning  before  gold  plating  was 
obtained  in  this  manner „ 
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A  frequently  troublesome  feature  of  (p,  # )  investigations  is  the 
gradual  accumulation  of  carbon  on  the  target,  as  a  result  of  cracking 
of  pump  oil  vapor  by  particle  bombardment.   To  reduce  this  difficulty, 
advantage  was  taken  of  an  8"in.  glass  cross  which  was  mounted  between 
the  target  diffusion  pump  and  the  target  chamber  as  shown  in  Fig.  5. 
A  special  liquid  air  cold  trap  was  constructed  and  installed  in  this 
cross  to  serve  a  dual  purpose.   First.,  it  functioned  as  an  ordinary 
cold  trap  with  a  large  surface  of  copper  vanes  suspended  from  the 
coolant  reservoir  to  provide  an  optically  dense  path  from  the  dif- 
fusion pump  thus  permitting  the  vacuum  to  be  routinely  maintained  close 
to  10  '  Torr.   The  second  feature  of  the  trap  was  a  1-in.  diameter 
copper  pipe  piercing  the  vanes  which  the  proton  beam  traversed.   At 
high  vacuum  the  proton  beam  functions  effectively  as  a  pump,  trans- 
porting gases  and  vapor  encountered  in  its  path.   When  the  beam  passes 
through  the  restricted  cold  pipe  for  6  in.  before  entering  the  target 
chamber,  it  was  reasoned  that  its  drifting  retinue  of  condensable  vapor 
contaminants  would  be  largely  trapped  by  the  cold  inner  wall.   The  re- 
sults from  use  of  the  trap  support  this  idea  since  only  minor  discolor- 
ration  was  observed  on  targets  compared  with  that  vnich  occurred  when 
the  trap  was  not  used,  and  no  increase  in  bombardment  energy  was  found 
necessary  as  targets  aged. 
4.   Shielding 

To  reduce  the  radiation  level  at  the  target  resulting  from  brems- 
strahlung  and  X-rays  caused  by  electrons  accelerated  toward  the  high 
voltage  accelerator  terminal 5  a  considerable  amount  of  lead  shielding 
was  erected.   A  2-in. -thick  lead  wall3  5  ft.  high  and  4  ft.  wide^  was 
placed  between  the  energy  control  slits  and  the  beam  collimator.   In 
addition,  two  lead  screens  1/4-in.  thick  and  each  5  ft.  by  4  ft.  were 
bolted  onto  a  caster  supported  steel  framework  constructed  for  this 
purpose.   The  mobility  of  the  screens  permitted  their  placement  to 
minimize  the  radiation  level  at  the  target  and  also  allowed  easy  re- 
moval for  access  to  the  target  area.   Although  the  above-described 

shielding  offered  some  reduction  in  natural  background  radiation, 

40 
(notably  the  1.46-Mev  K   gamma  ray)  it  was  not  feasible  to  screen 

specifically  against  background  without  either  Increasing  scattered 

target  radiation  at  the  detectors  or  utilizing  an  excessive  amount  of 

lead. 

10 


5.  Correlation  table 

A  correlation  table  shown  in  Fig.  6  was  built  for  this  experiment 
to  support  the  gamma-ray  detectors  and  to  permit  their  accurate  positioning 
with  respect  to  the  target.   A  steel  base  plate  mounted  on  casters  supported 
a  hollow  magnesium  column  upon  which  the  table  was  mounted.   Screw  jacks  in 
the  base  plate  were  extended  when  the  apparatus  was  in  place  to  hold  the 
table  in  position  and  to  allow  accurate  vertical  adjustment.   The  table  it- 
self  consisted  of  a  movable  graduated  aluminium  disc  which  was  mounted  con- 
centrically on  a  fixed  aluminum  support  plate.   Slides  mounted  on  the  disc 
and  on  the  support  plate  held  mounts  for  the  two  gamma-ray  detection  as- 
semblies.  The  mounts  were  positioned  in  the  slides  by  threaded  rods  which 
allowed  the  detectors  to  be  set  accurately  at  distances  from  1/2  in.  to 
6  in.  from  the  target  position.   The  detection  assembly  mounted  on  the  disc 
could  be  positioned  in  angle  with  respect  to  the  fixed  detector  in  12  min.  in- 
crements of  arc  by  means  of  an  adjusting  gear  meshing  with  gear  teeth  on 
the  periphery  of  the  disc.   The  entire  table  assembly  was  attached  to  the 
magnesium  support  column  by  a  hinge  so  that  it  could  be  tilted  into  either 
a  horizontal  or  vertical  plane s  permitting  measurements  to  be  taken  in 
either  the  horizontal  plane  containing  the  proton  beam  or  in  the  plane  per- 
pendicular to  the  beam.   The  first  named  arrangement  is  designated  herein 
as  geometry  A  and  the  latter  as  geometry  B. 

6.  Gamma-ray  detectors 

The  gamma-ray  detection  assemblies  were  supported  with  their  center- 
lines  8  3/4  in.  above  the  table.   They  each  consisted  of  a  3  in.  x  3  in. 
cylindrical  sodium  iodide  (thallium  activated)  crystal  integrally  mounted 
with  a  Dumont  6363  photomultiplier  tube  supplied  by  the  Harshaw  Chemical 
Company.   In  the  earlier  phases  of  this  work,  similar  2  in.  x  2  in.  crystals 
mounted  with  Dumont  6292  tubes  were  used.   The  resolution  of  the  3-in.  cry- 
stals, measured  after  the  signal  had  traversed  a  12-microsecond  delay  line 

137 
input  to  the  pulse  height  analyzer,  was  9%  for  the  0.662-Mev  Cs    gamma  ray. 

A  6AK5  cathode  follower  was  constructed  into  each  detector  mount  to  serve 

as  a  preamplifier. 

7.  Electronics 

The  electronic  equipment  and  circuitry  for  processing  signals  from 
the  detectors  were  planned  so  as  to  permit  considerable  flexibility  in  ex= 
perimental  applications.   A  block  diagram  of  the  basic  layout  is  shown  in 
Fig.  7. 

12 


o 


org 
5  in 

<  a 


cc  cc 
uj  o 
o  m 


T 

>iNva    a3Momod 
i 


3Q0H1V0 


^£&M&>* 


n 

I- 

o 

i 

z 

o 
cc 

cc 

hi 

t- 

cc 
u 

i- 
z 

< 

o 
o 

CO 

o 

z 
o 
cc 

\- 
o 
u 

-J 

UJ 


2 
< 

a: 
o 

< 


2*: 
o 
o 

CD 


UJ 

or 


14 


The  system  is  composed  of  an  RCL  128~channel  pulse  height  analyzer 
to  record  pulse  height  spectra,  a  slow-fast-sum  coincidence  analyzer,  and 
scalers  for  recording  various  outputs. 

It  was  found  that  requirements  for  the  detector  output  signal  by 
several  components ,  including  a  monitor  oscilliscope,  tended  to  overload 
the  preamplifiers  and  to  permit,  unacceptable  interaction  of  subsequent 
stages  with  each  other  through  the  coupling  of  their  inputs.   To  avoid  this 
problem  a  bank  of  six  cathode  followers  was  constructed.   A  bussing  arrange- 
ment made  it  possible  to  feed  the  grids  of  the  followers  from  either  one  of 
two  busses  or  from  an  individual  input,  so  that  up  to  six  pieces  of  equip- 
ment could  be  driven  bv  the  same  input  pulse  with  no  interaction  and  with 
practically  no  signal  loss.   Each  follower  was  provided  with  an  individual 
output  gain  control. 

The  output  of  either  detector  as  desired  was  also  fed  directly  into 
the  pulse  height  analyzer  through  a  simple  resistance  isolation  network. 
From  the  resistance  network  the  signal  passed  through  a  terminated  1350-ohm 
delay  line  to  allow  time  for  formation  of  gating  pulses  needed  for  coinci- 
dence spectra.   Several  methods  of  injecting  a  signal  into  the  analyzer 
were  tried  inasmuch  as  this  equipment  was  designed  to  accept  an  input  only 
from  a  photomultiplier  tube  connected  directly  into  the  preamplifier  as- 
sembly  provided.   This  was  not  an  acceptable  practice  because  there  existed 
no  convenient  means  of  signal  extraction  for  other  purposes,  and  because 
the  preamplifier  utiJ Lzed  a  high-voltage  supply  which  required  running  the 
photomultiplier  cathode  at  a  high  negative  potential  with  respect  to  the 
tube  shield  with  the  attendant  probability  of  tube  damage.   The  procedure 
adopted  as  most  successful  was  to  apply  the  signal  from  the  delay  line  to 
the  anode  pin  of  a  14-pin  tube  base  which  was  connected  directly  into  the 
analyzer  preamplifier.   tf  was  also  necessary  to  incorporate  a  gate  pulse- 
forming  network  into  the  analyzer  in  order  to  meet  the  gating  pulse  speci- 
fications required  by  this  equipment. 

Gamma-ray  spectra  were  obtained  by  gating  the  analyzer  with  the  dis- 
criminator output  pulses  from  the  linear  amplifier  fed  by  the  same  de- 
tector.  This  procedure  excluded  a  weak  precursor  pulse  which  was  trans- 
mitted by  the  delay  Lin<   aputs  by  not  gating  the  analyzer  on  until  a  few 
microseconds  before  arrival  of  the  delayed  pulse.   No  significant  loss  in 
detector  resolution  resulted  from  leaving  the  delay  line  in  the  input  for 
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ordinary  gamma-ray  spectra,  and  it  was  found  quite  convenient  to  be  able 
to  switch  from  these  singles  spectra  to  coincidence  spectra  without  change 
in  analyzer  calibration. 

For  the  observation  of  coincidence  spectra,  the  analyzer  was  gated  by 

the  slow-fast  coincidence  output*   The  basic  slow-fast  coincidence  system 

27 
assembled  for  this  experiment  has  been  described  elsewhere.     In  this 

system  the  output  from  an  Atomic  Instrument  Company  slow  coincidence  unit 

accepting  up  to  four  inputs  was  combined  in  a  diode  mixer  with  the  fast 

coincidence  unit  output  to  obtain  the  slow-fast  coincidence  output  signal. 

28 
The  fast  unit,  based  on  the  circuit  of  Draper  and  Fleischer,    detects  co= 

incidences  between  amplified  pulses  from  the  scintillation  detectors  with 

a  resolution  time  as  low  as  ten  ns.   For  pulses  spanning  a  large  amplitude 

range,  however,  more  dependable  operation  was  obtained  when  the  resolving 

time  was  increased  to  about  80  ns .   This  was  due  to  an  observed  amplitude 

dependence  of  the  shape  of  the  leading  edge  of  the  amplified  pulse,  to 

which  the  fast  unit  is  sensitive. 

Several  modifications  have  been  made  to  the  coincidence  system  which 

do  not  appear  in  the  previous  description.   The  diode  mixer  unit  for  com= 

bining  the  slow  and  fast  coincidence  outputs  has  been  modified  by  inserting 

a  Schmidt  trigger  between  the  adding  stage  and  the  output  stage  in  order  to 

clip  the  sum  pulse  and  improve  output-pulse  characteristics.   The  fast  unit 

has  been  modified  by  providing  a  variable  delay  line  in  the  input  in  order 

to  permit  accurate  compensation  for  differences  in  time  delays  between  the 

detectors  and  the  coincidence  unit.   Finally,  there  has  been  added  to  the 

system  a  provision  for  obtaining  sum-coincidence  spectra  as  first  described 

29 
by  Hoogenboom.     In  this  system  the  pulses  from  the  two  detectors  are 

carefully  matched  in  amplitude  calibration,  then  summed  in  a  pulse  adding 

circuit.   These  summed  pulses  are  analyzed  in  a  differential  discriminator 

set  to  accept  pulses  corresponding  to  the  excitation  energy  of  a  resonance 

level.   This  output  is  then  used  as  a  third  input  to  the  slow  coincidence 

unit  so  that  its  output  represents  coincidence  between  two  completely  ab~ 

sorbed  gamma  rays  comprising  a  double  cascade  from  the  resonance  level. 

This  signal  is  then  used  to  trigger  the  pulse  height  analyzer  gate.   The 

spectrum  obtained  from  pulse  height  analysis  of  one  detector  output  thus 

consists  only  of  full  energy  peaks  due  to  radiations  which  are  members  of 

a  double  cascade.   The  pulse  height  adder  which  has  been  constructed  for 

this  purpose  employs  the  circuit  used  by  Kantele  and  Fink  in  a  somewhat 

30 
different  type  of  spectrometer.    The  adder  output  was  analyzed  by  a  con- 
ventional linear  amplifier  and  discriminator. 
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The  scalers  used  to  record  single  channel  counts  had  dead  times  of  one 
microsecond  and  those  recording  coincidences  had  dead  times  of  approximately 
eight  microseconds.   Separate  cathode-f ollower  outputs  to  drive  scalers  were 
constructed  for  all  coincidence  stages  so  the  relatively  low  scaler  input 
impedances  would  not  degrade  coincidence  pulses.   The  number  of  scalers  pro- 
vided, which  exceeded  the  requirements  for  the  present  experimental  measure- 
ments, permitted  a  continuous  check  on  the  performance  of  the  coincidence 
system. 

A  double  pulser  with  associated  pulse  shaping  networks  and  an  os= 
cilloscope  formed  an  integral  part  of  the  equipment  and  were  used  for 

all  routine  checks  and  adjustments.   An  additional  slow  coincidence  unit 

31 
of  variable  resolving  time   ,  which  was  constructed  as  part  of  the  system, 

was  used  to  check  on  the  operation  of  the  Atomic  Instrument  Company  unit. 

Control  circuit  relays  were  provided  so  that  the  pulse  height  analyzer 

and  all  scalers  could  be  started  and  stopped  simultaneously  by  the  current 

integrator  on-off  switch. 
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III.   EXPERIMENTAL  PROCEDURE 

1.   Target  preparation 

The  target  material  used  in  this  experiment  was  isotopically  enriched 
silicon  in  the  dioxide  form  purchased  from  the  Stable  Isotopes  Division  of 

Oak  Ridge  National  Laboratory.   Relative  percentages  of  the  three  stable 

29         28  ^f) 

silicon  isotopes  were  91.85  Si   ,  7.85  Si   ,  and  0.30  Si   .   This  material 

was  evaporated  onto  prepared  target  backings  in  a  vacuum  evaporation  ap= 
paratus  under  a  chamber  pressure  of  10"   to  10   Torr . 

Beam  power  of  as  much  as  20  watts  was  delivered  to  a  target  spot  area 
of  about  6  square  millimeters  during  bombardment .   Dissipation  of  this  heat 
by  5-mil  tantalum  target  backings,  which  had  been  used  in  previous  work  at 
this  laboratory,  was  found  insufficient  to  prevent  local  overheating  of  the 
target  with  an  accompanying  rapid  loss  of  target  material.   This  fact,  and 
the  fluorine  contamination  of  the  surface  of  all  sheet  tantalum  examined, 
led  to  employment  of  the  heavier  target  backings  of  copper  described  in 
Chapter  II.   While  vacuum  plating  these  copper  backings  with  60,000  £  or 
more  of  gold  (99.999  percent  purity)  to  completely  stop  2-Mev  protons,  it 
was  found  that  the  plating  tended  to  crack  and  spall  off  when  the  copper 
backings  were  maintained  at  ambient  temperature  during  evaporation.   This 
effect  was  caused  chiefly  by  tensile  stresses  which  have  been  found  to  be 
as  much  as  850  Kg/cm   in  2000  £  gold  films  deposited  on  copper.32  This 
trouble  was  overcome  by  holding  the  copper  backings  at  400°C  during  eva- 
poration by  means  of  a  resistance  heater  which  was  constructed  for  this 
purpose.   The  heater  was  made  in  the  form  of  a  4=in.  cylindrical  section, 
to  which  the  target  backings  were  clamped.   The  cylindrical  heater  section 
was  suspended  in  the  vacuum  coating  chamber  so  that  its  axis  coincided 
with  the  evaporant  source.   Heating  the  backings  in  this  manner  prior  to 
evaporating  both  gold  and  Si02  also  served  to  remove  any  residual  volatile 
contaminants  leaving  a  clean  surface  to  receive  the  coating. 

Evaporating  boats  were  fabricated  from  5-mil  tungsten  ribbon  and,  when 
cleaned  by  preheating  before  use,  proved  to  be  satisfactory  for  the  eva- 
poration of  both  gold  and  Si02>   Temperature  of  the  Si0?  evaporant  was 
monitored  with  an  optical  pyrometer  to  maintain  an  indicated  evaporating 
temperature  of  2130  C;  at  this  temperature  4.0  mg  of  SiO,,  evaporated  in 
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about  two  minutes.   The  temperature  is  not  corrected  for  departure  of 
the  evaporant  from  an  ideal  black  body  inasmuch  as  reproducibility  and 
temperature  control  were  desired  rather  than  an  accurate  measurement  of 
true  surface  temperature.   During  the  first  five  to  ten  seconds  of  all 
evaporations  the  evaporant  was  intercepted  by  a  shutter,  since  the  first 
part  of  an  evaporant  to  boil  off  is  known  to  carry  a  higher  percentage 
of  any  contaminants  present .   Following  evaporation,  cooling  of  the 
heater  and  targets  was  accelerated  by  admitting  pure  helium  to  the  chamber 
through  a  liquid  air  trap,, 

In  practice,  2  gm  or  more  of  gold  were  evaporated  at  a  distance  of 
5  cm  from  the  target  backings  and  exposed  parts  of  target  assemblies  to 
produce  the  desired  plating  thickness.   4.1  mg  of  SiO~  were  evaporated  at 
the  same  distance  to  give  targets  approximately  5=kev  thick  to  2-Mev  pro- 
tons for  use  in  the  excitation  function  runs.   For  the  targets  used  to 
obtain  gamma-ray  spectra^  and  for  angular  correlation  runs  at  1470  kev, 
5.2  mg  of  SiO.  were  evaporated  to  obtain  targets  8-kev  thick  to  1470-kev 
protons. 

Considerable  care  was  used  in  an  attempt  to  reduce  the  contamination 

19  o 

by  F   .   Evaporation  boats  were  preheated  to  2200  C  before  use,  targets 

were  stored  in  a  vacuum  (until  this  was  found  not  to  be  of  value),  target 

chamber  pressure  was  maintained  at  10  '  Torr  during  runs,  and  possible 

sources  of  contamination  in  the  vacuum  systems  such  as  Teflon  were  a- 

19 
voided.   In  preliminary  bombardments  at  the  F   resonance  energy  of  1373  kev, 

the  gold-plated  copper  target  backings,  which  were  washed  in  either  pure 

ethanol  or  menthanol  and  rinsed  in  distilled  water  just  prior  to  use,  were 

19 
found  to  be  freer  from  F   than  backings  of  either  Ta  or  W  which  had  been 

cleaned  by  various  methods*   In  spite  of  these  precautions  all  targets 

19 

which  were  prepared  proved  to  be  contaminated  to  some  degree  by  F 

2.   Excitation  function  measurements 

29       30 
The  excitation  function  of  the  Si   (p,  K)P   reaction  has  been 

measured  by  recording  single  gamma-ray  yield  and  gamma-gamma  coincidence 
yield  at  accelerator  energy  settings  over  a  range  of  proton  energies  from 
1420  kev  to  2170  kev.   The  2  in.  x  2  in.  Nal(Tl)  scintillation  detectors 
were  used  for  this  purpose.   The  output  of  one  of  the  detectors  was  counted 
to  obtain  the  gamma-ray  yield  at  each  proton  energy,  and  the  number  of 
events  in  which  gamma  rays  were  detected  simultaneously  in  each  crystal  was 
recorded  by  the  use  of  the  slow-fast  coincidence  system  adjusted  to  a  re- 
solving time  of  25  ns . 
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The  two  detectors  were  positioned  in  geometry  B  at  angles  of  0  and 
180  relative  to  the  vertical  direction  and  at  a  distance  of  0.6  in.  from 
the  beam  spot.   Integral  discriminators  in  both  channels  were  set  so  that 
counts  would  be  recorded  for  gamma  rays  with  energies  exceeding  0.65  Mev, 

a  setting  chosen  to  exclude  the  0.51-Mev  annihilation  radiation  arising 

30 
from  the  positron  decay  of  P   ,  while  accepting  the  full  energy  peaks  of 

the  0.68  and  0.71  Mev  radiations  which  deexcite  the  first  two  excited 

30 
states  in  P   .   The  7-microampere  proton  beam  was  collected  by  the  cur- 
rent integrator  which  closed  a  beam  shutter  and  stopped  counters  after 
collection  of  84  microcoulombs  at  each  energy  step  of  approximately  3  kev, 
Background  counts  were  made  every  50  kev  with  the  beam  intercepted  by  a 
shutter  in  the  beam  tube  outside  of  the  shielded  target  area. 

An  auxiliary  experiment  was  conducted  to  determine  the  energy  cali- 
bration for  the  excitation  curves.   In  this  experiment  the  frequency  of 

"ft 
the  NMR  signal  was  measured  to  within  one  part  in  10   as  the  992-kev 

33 

aluminum  resonance    was  traversed  using  a  freshly  prepared  target  about 

10  kev  in  thickness.   Runs  with  mass  2  beams  were  also  made  to  check  this 

calibration  at  a  generator  energy  of  1984  kev.   Immediately  after  cali- 

29 
bration,  the  Si  '  resonance  energies  above  1970  kev  and  the  energy  of 

the  resonance  at  1506  kev  were  carefully  measured  in  a  similar  manner. 

Using  these  points  as  references  it  was  then  possible  to  calibrate  the 

entire  excitation  curve  in  terms  of  magnet  current,  which  was  determined 

to  a  precision  of  one  part  in  10   at  each  energy  step.   This  procedure 

gave  a  better  calibration  than  that  obtained  using  the  fluxmeter  logging 

dial  readings,  since  calibration  curves  for  converting  these  readings  to 

magnetic  field  could  not  be  read  with  adequate  precision. 

3.   Measurement  of  saturation  activity 

29       30 
In  order  to  definitely  identify  resonances  due  to  the  Si   (p,tf)P 

reaction,  measurements  were  made  at  the  resonances  above  1970  kev  to 

determine  the  half  life  of  decay  radiation  following  bombardment.   In 

29 
these  measurements  Si  '  targets  were  bombarded  for  ten  minutes,  then  re- 
moved to  a  small  cave  of  2-in.  lead  brick  where  counting  was  started 

within  60  sec  after  cessation  of  bombardment.   The  positron  annihilation 

30 
radiation  following  beta  decay  of  P   was  detected  by  a  3  in.  x  3  in. 

crystal.   Total  counts  above  the  noise  threshold  were  recorded  by  the 

pulse  height  analyzer  operated  in  a  gross  counting  mode.   In  this  mode 

the  channel  in  which  all  counts  are  recorded  is  automatically  shifted 

20 


ahead  after  a  preset  counting  time.   After  counting  for  3.1  sec  per 
channel  for  a  6.46-minute  total  counting  time,  the  target  was  removed 
from  the  cave  and  a  background  run  made  with  the  analyzer  in  the  sub- 
tract mode. 

Difficulty  which  was  experienced  in  obtaining  gamma-ray  spectra 

29 
recognizably  due  to  Si   at  the  upper  resonances  as  a  result  of  masking* 

radiation  from  the  F   (p,CC  "*T)0   reaction  caused  some  concern  that  the 

peaks  observed  in  the  excitation  curve  might  have  been  due  to  other  con~ 

30 
taminants.   Observance  of  the  2.5-min  half  life  of  P   following  bombard- 

ment  could  not  entirely  exclude  this  possibility  since  it  is  conceivable 

30 
for  P   to  be  formed  in  a  nonresonant  direct  capture  process.   For  this 

reason  it  was  desirable  to  compare  the  activity  of  the  target  observed 

after  bombardment  at  energies  corresponding  to  the  peaks  in  the  yield 

curves  to  the  activity  observed  after  bombardment  at  off-peak  energies. 

In  an  experiment  to  do  this>  a  proton  bombardment  time  of  ten  minutes  to 

obtain  94  percent  of  saturation  activity  was  followed  by  a  ten-minute 

counting  period  during  which  the  pulse  height  spectrum  was  obtained. 

The  target  was  then  allowed  to  decay  for  an  additional  ten  minutes  before 

bombardment  was  resumed  at  the  next  energy.   The  number  of  counts  observed 

in  the  0.51-Mev  annihilation  radiation  peak  then  served  as  a  measure  of 

30 
the  amount  of  P   formed  at  each  energy. 

4.   Gamma-ray  spectra 

All  gamma-ray  spectra  were  obtained  using  3  in.  x  3  in.  Nal(Tl) 
crystals  with  no  special  shielding  nor  collimation.   For  most  spectra 
the  pulse  height  analyzer  usable  range  extended  from  about  channel  7 
to  channel  127  corresponding  to  gamma-ray  energies  from  0.7  to  7.7  Mev. 
For  observation  of  annihilation  radiation  a  higher  gain  was  used  so  that 
these  channels  corresponded  to  energies  from  0.18  to  1.92  Mev.   The 
analysis  of  spectra  for  intensities  was  performed  with  a  digital  computer 
as  described  in  Appendix  2.   The  particular  method  finally  used  to  sub- 
tract background  radiation  from  spectra  varied  according  to  the  resonance 
being  considered.   When  possible,  backgrounds  taken  slightly  above  or  be- 
low resonances  were  subtracted  from  spectra.   This  was  not  very  successful 

19 
for  the  resonances  above  1950  kev  where  F   contamination  gave  rise  to  the 

F   (p,ajf)0   reaction.   This  reaction  is  strongly  resonant  at  2030  kev 

34 
with  a  width  at  half  maximum  of  120  kev   which  effectively  blankets  the 

higher  proton  energy  range  considered  in  the  present  work.   The  6-7  Mev 

21 


gamma  rays  occurring  in  this  reaction  are  present  in  relative  proportions 
which  vary  considerably  with  bombarding  energy  so  that  of f -resonance  back- 
ground  subtraction  is  not  accurate.   Bombardment  of  blank  targets  which 
had  been  processed  identically  to  SiO.  targets  but  using  an  empty  eva- 
poration boat  gave  a  background  measurement  which  was  also  rather  unre- 
liable because  of  individual  variations  among  targets.   A  method  of  sub= 
tracting  the  background  measured  with  the  accelerator  off  was  also  used, 
in  which  case  background  radiation  of  specific  energies  due  to  contaminants 
was  subtracted  out  in  the  spectrum  analysis. 

In  order  to  circumvent  radiation  due  to  target  contaminants  efforts 
were  exerted  to  obtain  sun-coincidence  spectra  in  the  energy  region  above 
1950  kev  as  a  means  of  getting  at  decay  schemes.   The  sum-coincidence 
spectra  were,  in  general,  measured  with  the  slow-fast  coincidence  require- 
ment  also  imposed,  and  with  a  coincidence  resolving  time  of  110  ns.   These 
runs  were  usually  made  in  geometry  B  with  counters  positioned  1.6  in.  from 
the  beam  spot  at  a  90  angle  relative  to  each  other.   A  5/16-in.  lead  sheet 
between  the  counters  served  to  minimize  coincidences  due  to  scattering  or 
pair  production  in  one  crystal  with  the  lost  radiation  being  captured  in  the 
other  crystal.   The  differential  discriminator  following  the  pulse  adder  was 
set  with  its  lower  level  at  the  resonance  energy  and  was  adjusted  to  pass 
pulses  in  a  band  about  0.6  Mev  in  width.   Settings  were  made  by  adjusting 
the  discriminator  so  that  the  observed  count  rate  of  accepted  pulses  was 
equal  to  half  the  repetition  rate  of  a  test  pulse.   The  proper  test  pulse 
amplitude  was  obtained  by  setting  the  pulser  gain  so  that  counts  fell  into 
the  pulse  height  analyzer  channel  corresponding  to  the  desired  energy. 
Settings  of  integral  discriminators  in  the  two  singles  channels  were  adjusted 
to  0.62  Mev  in  a  similiar  manner.   The  pulse  height  analyzer  was  calibrated 

with  Cs    and  Co   sources  and  with  the  5.31-Mev  gamma  ray  which  dominates 

29       30 
the  Si   (p,  }f)P   spectrum  at  the  414-kev  resonance.   It  was  found  that  the 

Dumont  6363  photomultiplier  tubes  experienced  marked  gain  shifts  as  a  function 

of  count  rate  amounting  to  several  percent,  an  effect  which  has  been  attri- 

35 
buted  to  changes  in  the  dynode  surfaces   .   Attempts  to  allow  this  shift 

to  stabilize  before  commencing  counts  were  only  partially  effective.   The 

overall  result  of  this  phenomenon  was  to  vitiate  the  results  of  a  large 

number  of  sum-coincidence  spectra. 
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5.   Measurements  at  1470  kev 

30 
Since  the  4.50-Mev  level  in  P   has  been  shown  to  be  most  strongly 

excited  at  the  1470-kev  resonance,  this  proton  energy  was  a  logical  choice 
for  angular  correlation  measurements  which  could  lead  to  determination  of 
the  spin-parity  for  the  level.   The  resonance  is  rather  weak  so  it  was  de- 
cided to  take  five  sets  of  data  to  improve  the  counting  statistics  and  to 
permit  a  check  of  the  consistency  of  results.   Each  set  of  data  consisted 
of  five  spectra  taken  with  the  movable  detector  set  to  angles  of  0  ,  30  , 
45  ,  60  ,  and  90  relative  to  the  proton  beam.   In  geometry  A  with  the 
3  in.  x  3  in.  movable  detector  positioned  with  its  face  4.125  in.  from 
the  target,  runs  were  made  using  a  7 -microampere  beam  for  counting  times 
of  approximately  one  hour  at  each  angular  setting.   The  same  target  was 
used  for  each  set  of  runs  and  the  order  of  the  angular  settings  was  varied 
in  different  sets  so  as  to  minimize  differences  in  the  gamma-ray  spectra 
due  to  any  contaminant  buildup  or  target  erosion.   The  fixed  3  in.  x  3  in. 
detector  was  positioned  at  270  relative  to  the  beam  and  4.125  in.  from 
the  target  for  use  as  a  monitor  counter.   Each  run  was  terminated  when 
30,000  monitor  counts  had  been  recorded.   A  2  in.  x  2  in.  detector  was 
also  used  as  a  secondary  monitor  counter.   The  integral  discriminator  in 
each  monitor  channel  was  adjusted  to  a  valley  in  the  gamma-ray  spectrum 
at  4.75  Mev  so  that  the  effect  of  gain  shifts  in  these  channels  would  be 
minimized.   Before  and  after  each  run  the  discriminator  levels  were  checked 
and  adjusted  with  the  aid  of  a  Co   source  which  was  precisely  and  repro- 
ducibly  located  relative  to  the  counters.   For  this  check  the  gain  in  the 
monitor  channels  was  increased  by  a  factor  of  4  by  means  of  the  gain  mul- 
tiplier switch  so  that  the  discriminator  setting  was  equivalent  to  1.14  Mev. 
At  this  setting,  near  the  1.17-Mev  peak  in  Co   ,  the  count  rate  observed 
with  the  Co   source  was  quite  sensitive  to  the  discriminator  setting  so 
that  accurate  adjustment  was  possible. 

The  unfolding  of  gamma-ray  spectra  and  the  determination  of  angular 
distribution  coefficients  was  programed  and  analyzed  on  a  CDC-1604  computer 
as  described  in  Appendix  2.   The  radiations  which  have  been  used  as  spectro- 
meter calibration  curves  in  this  analysis  are  listed  in  Table  1.   These 
spectra  were  measured  in  experimental  geometry  As  backgrounds  were  sub- 
tracted, and  any  remaining  contaminating  gamma  rays  were  subtracted  off. 
Each  of  these  monoenergetic  spectrograms  was  plotted,  hand  smoothed,  and 
redivided  into  525  channels  to  reduce  any  interpolation  errors  in  the 
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TABLE  1 
Gamma  Rays  Used  for  Spectrometer  Calibration 


E v(Mev)  Radiation  Source        E  (kev) 

t  P 


0.662 

1.63 

2.37 

2.75 

3.51 

4.43 

5.31 

6.135 

8.06 


Cs137 

- 

„  23,    r.„   20 
Na   (p,a  0)Ne 

1252 

12       13 

460 

«  24 

Na 

- 

_12,   -.v„13 

c   (p,  dh 

1698 

1 1       12 

675 

Si29(p5(T)P3() 

414 

^19,    ^^16 
F  (psa(T)0 

1374 

13       14 

c    (p,  *r)N 

554 

24 


subsequent  analysis.   Energy  calibration  of  each  spectrum  was  made  by 
equating  distance  in  channels  between  the  pair  peak  and  the  photopeak 
to  1.022  Mev.   For  the  three  lowest  energy  calibration  curves,  for  which 
pair  peaks  are  non-existent  or  poorly  defined,  an  equivalent  calibration 
was  computed  based  on  the  pulse  height  analyzer  calibration  as  determined 
in  an  auxiliary  experiment . 

All  spectra  were  read  out  of  the  pulse  height  analyzer  memory  photo- 
graphically.  These  data  were  then  transcribed  to  punched  cards  and  printed 
listings  of  the  cards  were  carefully  proofread  against  the  photographs. 

Gain  drifts  did  not  prove  troublesome  at  the  lower  count  rates  ob- 
served during  correlation  runs  at  1470  kev.   To  compensate  for  differences 
in  calibration  of  the  pulse  height  analyzer  between  successive  runs,  which 
were  of  the  order  of  one  percent,  each  spectrum  was  self -calibrated  by  the 
peaks  at  7o00  Mev  and  at  1.45  Mev. 
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IV  EXPERIMENTAL  RESULTS 

1.   Excitation  Function 

29       30 
The  excitation  function  for  Si   (p,}f)P   as  measured  by  gamma-ray 

yield  is  shown  in  Fig.  8,  and  as  measured  by  gamma-gamma  coincidence 
yield  in  Fig.  9.   The  energies  at  which  resonances  are  observed  are  listed 
in  Table  2,  together  with  some  of  the  results  which  have  been  obtained  at 
proton  energies  below  1900  kev  by  other  investigators. 

In  arriving  at  the  uncertainty  of  0.207o  kev  for  the  resonance  energies, 
consideration  has  been  given  to  errors  arising  from  instability  and  inhomo- 
geneity  of  the  field  of  the  analyzing  magnet,  the  proton  beam  spread,  pos- 
sible surface  contaminants,  calibration  drift,  and  the  0.5-kev  estimated 

33 
error  for  the  992-kev  aluminum  calibration  point. 

Particular  attention  was  given  to  the  resonances  above  1900  kev  which 

have  not  previously  been  reported.   The  results  of  determinations  of  the 

half  life  of  the  post -bombardment  activity  of  the  targets  at  these  energies 

are  given  in  Table  3.   The  measurements  are  seen  to  agree  with  values  for 

30 
the  half  life  of  P   reported  in  the  literature,  which  range  from  2.51  + 

24  29 

.02  min  to  2.62  +  .02  min,    so  that  radiative  capture  by  Si  '  is  clearly 

demonstrated.   Table  4  shows  the  intensity  of  positron  annihilation 

radiation  observed  after  saturation  at  selected  bombardment  energies. 

These  intensity  values  show  distinct  maxima  at  1975,  2033,  2075,  and  2117 

30 
kev  and  show  that  P   is  formed  in  a  resonant  reaction  rather  than  by  a 

non-resonant  process.   The  peaks  in  the  gamma-ray  yield  curves  at  these 

energies  can  therefore  be  positively  ascribed  to  resonant  proton  reaction 

29 
with  Si   . 

The  yield  of  0. 511-Mev  radiation  after  bombardment  at  energies  between 

resonances  is  probably  due  to  a  combination  of  off-resonant  reaction  with 

29  17 

Si   and  to  the  66-sec  positron  activity  of  F   following  the  non-resonant 

0   (p,  If )F   reaction.   This  hypothesis  is  supported  by  half  life  measure- 
ments of  post-bombardment  activity  of  the  targets  at  two  inter-resonance 
points  at  1955  and  2007  kev»   At  those  energies  "half -lives"  of  the  order 
of  1.6  +  .5  min  are  observed  with  the  count  rate  2  min  after  bombardment 
at  1955  kev  a  factor  of  ten  lower  than  the  count  rate  2  min  after  bombard- 
ment at  1975  kev. 
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TABLE  2 

29       30 
Resonances  Observed  in  the  Si   (p,  if)?       Reaction 

E  from  1450  to  2170.   Values  in  kev. 
P 

This  work   Green     Seagondollar     Tsytko     Wild     Remarks 


1470 

1479 

1479 

1470 

1506 

1515 

1512 

1505 

1637 

1648 

1643 

1635 

1638 

(1646) 
1667 

1671 

(1653) 

(1647) 
1663 

(165 
(1665 
(1670 

1688 

1692 

1693 

1680 

1685 

(1726) 

1752 

1752 

1752 

1749 

1775 

1777 

1775 

1772 

(1803) 

1811 

1808 

1860 

1857 

1852 

(1933) 

1975 

(2007) 

2033 

2075 

2117 

(2159) 

Si 


28 


weak 


weak 


Si 


30 


not  identified 


+0 . 20% 


+2 


+2 


+3 


Errors 


29 


TABLE  3 

30 
Half  Life  of  Post -bombardment  Decay  of  P 


Bombardment  Energy  (kev)  Half  life  (min) 

1975  2.42  +  .09 

2033  2.45  +  .18 

2075  2.51  +  .06 

2117  2.47  +  .15 
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TABLE  4 
Intensity  of  Post-bombardment  Annihilation  Radiation 


Relative 
Intensity 

19,3 

1.0 

4.1 

1.0 

7.9 

1.8 
12.4 

1.3 

8.1 

2.3 


Bombarding 
Energy  (kev) 

Excitation  «_ 
Energy  in  P   (Mev) 

1506 

7.04 

1950 

- 

1975 

7.49 

1990 

= 

2033 

7.55 

2060 

- 

2075 

7.59 

2100 

- 

2117 

7.63 

2140 

_ 

* 
Based  on  Q  =  5.582,  reference  24. 
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2.   Gamma-ray  Spectra  at  1975,  2033,  2075  and  2117  ksv  Resonances 

Gamma-ray  spectra  were  obtained  at  the  four  resonances  at  proton 

energies  of  1975,  2033,  2075  and  2117  kev.   Masking  radiation  from  the 

19 
resonant  reaction  with  the  contaminant  F   dominates  the  upper  energy 

region  of  these  spectra  so  that  a  detailed  analysis  of  the  relatively 

29       30 
weaker  radiation  from  the  Si   (p, ft  )P   reaction  was  not  possible. 

However,  the  gamma-ray  energies  and  intensities  which  could  be  deter- 
mined with  reasonable  certainty  from  stripping  and  least  squares  analysis 

19 

for  intensities  are  given  in  Table  5.   The  6=7  Mev  gamma  rays  due  to  F   , 

which  are  omitted  in  Table  5,  were  accounted  for  by  inclusion  in  the 
least  squares  intensity  analysis,  since  other  background  subtraction 
techniques  did  not  work  well  at  these  energies  as  has  been  discussed 
in  Chapter  III.   Also  included  in  Table  5  are  energies  of  sum-coinci- 
dence spectral  peaks  which  were  close  to  energies  observed  in  the 
singles  spectra.   The  good  spectral  resolution  and  discrimination 
against  background  radiation  which  are  obtainable  with  the  sum-coinci- 
dence method  proved  to  be  of  only  nominal  assistance  in  analyzing  the 
gamma-ray  spectra.   This  was  due  in  part  to  the  relative  weakness  of 
radiative  decay  of  these  resonance  levels  in  comparison  with  inelastic 
scattering,  and  in  part  due  to  photomultiplier  gain  drifts  which  blurred 
peaks  in  the  sum-coincidence  spectra.   However,  it  is  apparent  from  the 
large  number  of  low  intensity  peaks  observed  in  the  sum-coincidence 
spectra  at  the  higher-energy  resonances  that  decay  schemes  from  these 

levels  are  considerably  more  complex  than  decay  from  the  1752-kev  re- 

23 
sonance  level,    where  sum-coincidence  spectra  were  recorded  for  the 

purpose  of  comparison. 

The  strong  1.28-Mev  gamma  ray  observed  at  1975,  2033,  and  2075  kev 

30 
is  interpreted  as  deexcitation  of  the  P   compound  nucleus  by  proton 

29  25 

reemission  to  the  first  excited  state  of  Si  '  at  1.277  Mev.    The  high 

intensity  of  this  gamma  ray  is  convincing  evidence  that  this  inelastic 

scattering  process  predominates  over  radiative  capture  and  that  the 

128-Mev  radiation  is  primarily  responsible  for  the  upper-energy  peaks 

in  the  gamma-ray  yield  curve.   Rather  curiously,  however,  this  gamma 

ray  does  not  appear  in  the  spectrum  at  the  2117-kev  resonance. 

The  resonance  at  2075  kev  lies  near  the  Si   (Pj^P  '  resonance 

25 
reported  at  2090  kev  which  has  a  half  width  of  16  kev.   The  gamma-ray 
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TABLE  5 

Gamma-ray  Intensities  at  Resonance 
Energies  of  1975,  2033,  2075,  and  2117  kev 


E   (kev) 
P 

Singles  Spectra 
E  „  (Mev) 

Sum-Coincidence 
Ey(Mev) 

1975 

4.52 

4.58 

4.19 

4.07 

3.48 
3.33 

3.40 

3.00 

2.99 

2.44 

2.50 

1.28 

- 

2033 

5.10 

- 

1.28 

- 

2075 

7.59 

- 

4.70 

4.85 

2.98 

2.95 

2.26 

2.23 

1.36 

- 

1.28 

- 

2117 

5.10 

- 

2.55 

= 

Intensity  (quanta/A- C) 


Error  in  singles  spectra  E^  estimated  to  be  1% 
Error  in  sum-coincidence  E*  estimated  to  be  2% 
Error  in  intensities  estimated  to  be  30% 


5 
2 

2 
2 

3 

3 

117 

8.6 

187.3 

14 

60 

55 

55 

303 

1360 

2 

3 
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TABLE  6 

Relative  Gamma-ray  Intensities  at  the  1470-kev  Resonance 

_ 

E  v»  (Mev)    Transition         Observed  Decay  Scheme   Ref  23 

Intensity  Intensity 

7.00       R->0               25.4  25.4  39.6 

6.29       R-^>  2                9.6  9.6  6.3 

4.50       19-^0              14.6  14.3  13.2 

4.14       13-*0              10.7  10.7  11.1 

3.98       R-»-9               10.5  10.8  7.3 

3.45       not  accounted  for      2.0  =  4.0 

3.17       R-=»ll  1.5  2.0 

3.17       ll-*.l  1.5  2.0 

3.07       R-^.12  17.6  11.7 

3.02       9-*0  1.2  1.6 

2.86       R-»-13              13.7  13.7  13.0 

2.69       13-*-3               2.8  2.8 

2.50       R-p-19                QO  0  21.6  20.0 

JO  .  o 

2.50       12-*- 3  17.6  11.7 

2.34       9-»-l               17.2  16.9  13.4 

1.45       19-**  9            not  analyzed"1"    7.3  6.8 

1.45       3-*-0                "  20.3  11.7 


.71       2->0 
.68       l->0 


21.8 


An  intensity  of  approximately  34  is  indicated  for  the  1.45  Mev  radiation 
by  observations  at  90°  with  .6"  crystal  face  to  target  distance. 

Relative  intensities  shown  are  also  approximate  absolute  intensities 
in  quanta  per  microcoulomb. 


34 


3 


CD 

o 

CO 

o 
en 

<u 
Of! 

> 
CD 
J»5 

0 
O 


CO 

u 


14-1 

CD 
O 

o 

g 

o 

t4 
4J 

CO 
.-1 
CD 
U 

u 

o 
o 

t-l 

CO 
1-1 

3 
60 

d 

< 


3 
O 


CO 

CD 

X! 
60 


CO 


4-1 
CO 
CD 

XI 
00 

•l-l 


CN 


CN 


co 
o 


lO 

o 


o 


on 
o 


m 
o 


O 


O 


CN 

o 


in 
o 


o 


ON 


CN       •* 
O      O 


i-l      CO 
CO      o 


o 


CO 

o 


o 


CN 

o 


o     o 


o 


in     cm 

r-l         O 


CO 


00 

o 


CO 


CN 

CO 


o 
m 


CN 


CO 

o 


CN 


CM 
O 


m 
o 


CN 

o 


o 


m 
o 


CN 
CM 


vO 


O 


CO 

O 


CM        CO 

O      O 


CM 

r^ 

ON 

ON 

m 

co 

Ov 

<f 

<t 

<T 

o 

O 

CO 

o 

o 

o 

CO 


on 

CO 


CN 
CN 


CM 


> 

CD 

X 

s_/ 

o 

o\ 

o 

•So 

o 

CM 

m 

M 

r-*- 

vO 

<f 

00 

ON 


m 


r^  cm 

vO 

on 

O 

<f 

O  O 

00 

sO 

in 

CO 

u 

o 

u 

CD 
4-1 
CD 
"0 


C 
o 
•1-1 

4_> 

3 

r-l 
O 

co 
CD 
l-i 

u 

CO 
i-l 

3 
00 
C 
CO 


a 

•H 

M-l 


T3 
CD 
4-1 
CJ 
CD 
U 

u 

o 
o 

CD 
J-i 
CO 


O 

Xi 
to 

CO 

4-1 

c 

CD 

•I-l 
o 
•I-l 

14-1 
CM 
CD 

o 
o 


COCOCOCOCOCMCMCMCN 


35 


spectrum  at  this  resonance  is  characterized  by  a  4.76-Mev  radiation 
to  the  ground  state  of  P        It  seems  probable,  therefore,  that 

the  4.70-Mev  gamma  ray  seen  in  the  spectrum  at  2075  kev  is  due  to  the 

28 
Si   reaction.   The  observed  intensity  is  in  agreement  with  this  ex- 
planation.  Although  the  presence  of  2.94-Mev  and  2.24-Mev  gamma  rays 

30 
in  this  spectrum  suggest  a  cascade  through  level  8  at  2.94-Mev  in  P   , 

22 
the  1.45-Mev  gamma  ray  known  to  occur  in  this  cascade   is  not  seen 

and  sum-coincidence  work  has  not  been  able  to  confirm  this  possibility. 

29 
It  is  possible,  however,  that  the  gamma  rays  due  to  Si   reactions  could 

be  responsible  for  the  anisotropy  of  the  radiation  from  the  4.76-Mev 

29  37 

level  in  P   which  has  been  observed  by  Newton,    since  anisotropy  is 

not  expected  to  result  from  the  spin  assignment  to  this  level  of  1/2+ 

tJ  38 

based  on  proton  scattering  experiments. 

The  remaining  weaker  radiations  seen  at  the  four  higher  resonances 

29 
are  most  probably  due  to  radiative  capture  by  Si   ,  but  the  rather  large 

errors,  due  chiefly  to  masking  contaminant  radiation,  do  not  permit  more 

than  tentative  postulation  of  decay  modes.   Gamma  rays  observed  at  the 

1975-kev  resonance  suggest  cascades  through  level  19  at  4.50  Mev  and 

through  level  14  at  4.18  Mev.   No  conclusions  can  be  drawn  for  the 

2033-kev  resonance.   At  the  2075-kev  resonance  the  7.59-Mev  radiation 

is  readily  recognized  as  a  ground  state  transition  but  the  other  gamma 

rays  are  not  identified.   Finally,  at  the  2117-kev  resonance,  the  two 

radiations  observed  nicely  fit  a  decay  through  level  5  at  2.54  Mev  which 

22 
is  known  to  decay  95%  to  the  ground  state. 

3.   1470-kev  Resonance  Decay  and  Intensities 

The  gamma-ray  spectrum  observed  at  the  1470-kev  resonance  with  a 
crystal  face  to  target  distance  of  0.6  in.  and  at  an  observation  angle 
of  90   is  shown  in  Fig.  10.   This  spectrum  has  been  least  squares  fitted 
by  gamma  rays  with  energies  required  by  the  decay  scheme  proposed  in  ref- 
erence 23.   Reanalyses  with  alternate  energy  choices  in  each  instance  gave 
a  poorer  fit  in  reassurance  of  the  correctness  of  the  decay  modes  proposed, 
This  decay  scheme,  with  intensities  as  found  in  the  present  analysis,  is 
shown  in  Fig.  11.   (Heavy  dots  indicate  cascades  confirmed  in  reference  23 
by  coincidence  spectroscopy.)   The  branching  ratio  shown  for  the  4.50-Mev 
level  is  that  given  in  reference  22.   The  assumed  correctness  of  this 
ratio  then  requires  branching  ratios  of  93:7  for  level  9  and  80:20  for 
level  13  in  order  to  fit  the  intensities  of  gamma  rays  which  have  been 
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DECAY     SCHEME      AT      1470      KEV      RESONANCE 
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observed.   These  ratios  are  in  reasonably  good  agreement  with  ratios  of 

22  23 

89:11  reported  for  level  9   and  88:12  reported  for  level  13 

Intensities  of  gamma  rays  observed  in  the  angular  correlation  runs 

are  shown  in  Table  6  along  with  intensities  observed  in  reference  23  at 

55  .   The  intensities  found  in  the  correlation  runs  have  been  averaged 

over  4jt  steradians  by  weighting  the  mean  values  for  each  observation  angle 

by  the  sine  of  that  angle .   Because  the  1.46-Mev  background  radiation  due 

40 
to  K   was  found  to  be  quite  strong  relative  to  the  1.45-Mev 

gamma  ray  observed  in  the  correlation  runs,  no  attempt  was  made  to  analyze 
the  gamma  radiation  of  this  or  lower  energy.   Although  the  background 
radiation  subtracts  out,  the  statistics  of  the  remaining  counts  do  not 
allow  extraction  of  significant  information.   The  observed  intensities 
are  seen  to  agree  reasonably  well  with  those  of  reference  23  except  for 
the  7.00=-Mev  radiation  which  is  found  to  be  considerably  weaker  in  the 
present  work.   Since  the  experimental  intensities  are  based  upon  reported 
values  of  intrinsic  peak  efficiencies   rather  than  an  absolute  spectro- 
meter calibration,  this  discrepancy  is  not  considered  significant. 
4.   Angular  Correlations  at  the  1470-kev  Resonance 


A.   General 

The  correlation  function  coefficients  which  have  been  found  in  this 
experiment  are  given  in  Table  7.   The  values  given  are  weighted  means  of 
the  values  found  in  each  of  the  five  sets  of  data  taken.   The  chi-square 
values  for  each  correlation  coefficient  have  been  calculated  and  show  the 
different  sets  to  be  consistent,  under  the  acceptance  requirement  that  the 
cumulative  chi-square  probability  be  equal  to  or  less  than  0.9.   Combining 
intensities  from  each  data  set,  then  analyzing  for  the  coefficients  gives 
values  practically  identical  to  the  mean  values  quoted.   This  agreement 
to  some  degree  serves  as  a  check  that  computational  errors  have  not  crept 
in.   Each  experimentally  observed  spectrum  has  also  been  plotted  together 
with  the  corresponding  fitted  spectrum,  then  visually  examined  to  determine 
that  a  reasonable  fit  was  obtained  in  each  case.   Fig.  12  shows  a  typical 
example  of  the  fit  obtained  for  the  correlation  runs. 

The  intensity  analysis  was  performed  both  subtracting  background 
radiation  analytically  and  subtracting  of f -resonance  background,  with 
essentially  the  same  results  being  obtained  in  each  case.   The  standard 
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deviations  in  the  mean  coefficient  values  have  been  calculated  from 

i=l 
and  have  been  compared  with  the  values  given  by 


1     " 
n(n-l)  Z  C 


n  « 

x.  -  x) 


which  agree  closely, 


The  most  obvious  feature  of  the  data  is  the  anisotropy  of  the 
4.50~Mev  gamma  ray  which  is  incompatible  with  a  spin  of  zero  for  this 
level.   It  is  also  noteworthy  that  the  4.50-Mev  gamma  ray  is  the  only 
radiation  to  show  a  statistically  significant  P. (cos  9)  term  (notation 
is  defined  in  Appendix  1) .   Fig.  13  shows  the  correlation  coefficient 
values  for  this  gamma  ray  obtained  from  each  set  of  data.   In  this 
figure  M  denotes  mean  values  and  C  denotes  values  obtained  by  com= 
bining  intensities  from  each  data  set  before  computing  coefficients. 
Fig.  14  shows  the  fit  to  the  Legendre  series  both  with  and  without  P, 
terms.   Although  the  existence  of  the  P  term  is  not  incontrovertible 
it  permits  a  markedly  better  fit  to  the  data  and  gives  a  smaller  value 

of  the  parameter  £  which  serves  as  a  guide  to  the  appropriate  number  of 

41 
terms  which  should  appear  in  the  Legendre  series. 

In  order  to  deduce  spin-parity  values  it  is  necessary  to  examine 

the  possible  particle  and  radiation  mixing  ratios  for  this  reaction. 

The  notation  used  herein  to  refer  to  these  mixing  ratios  is  as  follows: 

3-,  -  Proton-partial-wave  mixing  parameter 

g_  =  The  multipole  mixing  ratio  of  an  intermediate 
unobserved  gamma  ray 

r_  =  The  multipole  mixing  ratio  of  an  observed 
gamma  ray 

t,    =  Channel-spin  mixing  parameter 

Appendix  1  defines  these  parameters  more  fully  and  describes  how  they 

enter  into  the  expression  for  the  theoretical  correlation  function.   For 

29       30       29        30 
the  Si   (p5tf)P   and  Si   (p,x(f)P   (intermediate  unobserved)  reactions, 

41 


.20 


.10 


A2  /Aq 


0 


9      -r- 


6 


h 


-.10 


-.20 


A4 

/A 

0 

T 

( 

) 

j 

IT 

6 

L  1  I 

T 

6 

RUN  12     3     4     5       M       C 


P4    HIGHEST    TERM        6  =    .74 


2     3     4     5       M       C 


RUN 


.20 


.10 


2     3     4     5      M       C 

A2/A0 


0 


6 


0 


H 


Pz    HIGHEST    TERM      €  =  1.30 

FIGURE     13.     DISTRIBUTION     COEFFICIENTS      FOR    THE 
4.50    MEV     GAMMA      RAY 


42 


1.20    * 


I   +    .09P2    (COS  9)     -     .05P4    (COS  6) 
I    +  .09P2    (COS  90°)-    .05P4    (COS  90°) 


1.10 


z 

£   1.00    - 


1.20    " 


€  *   0.74 


o  o  o 

30         45  60 


e 


90" 


1.10 


g  1.00    - 

UJ 


I     +     .07P2   (COS0) 


6  ■  1.30 


30°         45 

e 


US 


I    +     .07  P2    (COS   90) 


\ 


60' 


90 


FIGURE        14.  FITS   OF    COMBINED     EXPERIMENTAL     DATA      TO 

LEGENDRE      SERIES    -    4.50-MEV     GAMMA      RAY 
(UNCORRECTED      FOR    FINITE     SOLtD     ANGLE) 


43 


the  correlation  coefficients  A0/A~  and  A. /A„  are  describable  as  functions 

^  0      4  0 

of  a  maximum  of  three  mixing  ratios.   Level  curves,  or  contours,  of  con- 
stant A./An  values  may  be  plotted  against  any  two  of  the  mixing  ratios. 
Actually  the  arctangents  of  the  mixing  ratios  are  used  as  variables  for 
convenience  since  they  remain  finite.   If  a  third  mixing  ratio  appears, 
changing  its  value  changes  these  contour  plots  only  to  the  extent  of 
changing  the  A./A0  values  for  all  contours  by  a  multiplicative  constant. 
Contour  plots  for  this  analysis  have  been  programmed  and  plotted  as 
described  in  Appendix  2. 

In  the  following  discussion  it  will  be  assumed  that  only  El  radiation 
or  E2/M1  mixtures  are  likely  to  be  observed  for  transitions  from  levels 

of  several  Mev.   This  simplifying  assumption  is  based  upon  transition 

42 
probabilities  estimated  from  the  independent  particle  model    and  is  com- 
monly made  in  analysis  of  gamma  radiation  following  proton  capture  in 
light  and  intermediate  nuclei..   Although  suppression  or  inhibition  of  di- 

pole  radiation  is  possible  because  of  the  isotopic  spin  selection  rule, 

a  43 

&  T  =  +  1  for  self  conjugate  nuclei,    there  are  an  adequate  number  of 

30 
lower  levels  of  both  T  =  0  and  T  =  1  in  P   so  that  the  resonance  level 

can  decay  by  either  El  or  E2/M1  radiation  without  violating  this  rule. 

In  analyzing  the  experimental  correlation  coefficient  values,  tables 

from  reference  23  have  been  quite  helpful.   This  reference  tabulates  the 

maximum  and  minimum  values  which  may  be  assumed  by  A~/An  and  A, /An  for  a 

variety  of  postulated  spin  sequences. 


B.   Resonance  level  spin-parity 

The  observed  A  /An  terms  of  about  -0.40  for  the  7.00,  6.29,  and 
3.98  Mev  gamma  rays,  which  all  proceed  to  known  1+  levels,  are  found  to 
be  most  useful  in  determining  the  spin-parity  of  the  resonance  level. 
Among  possible  values  through  4  for  this  level,  it  is  found  that  only 
1+,  2+,  and  2-  can  lead  to  the  A~/An  coefficient  measured  for  the  7.00-Mev 
ground  state  transition.   The  same  conclusion  is  supported  by  the  A~/An 
values  observed  for  the  6.29  and  3.98  Mev  radiations. 

If  one  inspects  the  contour  plot,  Fig.  15,  for  the  1+  — >  1+  transition 
it  is  seen  that  a  d-s  proton-angular -momentum  mixing  ratio  of  over  0.90 
is  required  to  produce  the  observed  A„/An  value  for  the  7.00-Mev  gamma 

ray.   Since  the  relative  penetrability  of  d  waves  to  s  waves  is  about 

-2 
3  x  10   ,  such  a  mixture  is  not  a  likely  choice.   The  A,/An  term  observed 
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for  the  4„50-Mev  radiation  supports  this  conculsion  since  A . /A„  terms 
do  not  result  from  resonance  spins  smaller  than  2.   Finally,  the  value 
of  A_/An  observed  for  the  2.34-Mev  radiation  will  be  shown  to  be  incom- 
patible with  a  1+  resonance  level  spin.   This  radiation  proceeds  to  the 
known  0+  level  1,  at  0.68  Mev,  from  the  known  1+  level  9, at  3.02  Mev. 
Level  9  is  fed  more  than  507o  from  the  resonance  level  and,  secondarily, 
by  the  intermediate  1.48-Mev  radiation  from  level  19.   The  A„/An  value 
observed  for  this  primary  3.98-Mev  feeder  restricts  the  ranges  of  ^ 
and  S9  for  the  R— *9— *  1  transition.   When  these  ranges  are  delineated 
on  a  contour  plot  for  a  1+  -^-1+  -*0+  transition,  it  is  found  that  A9/An 

for  the  second  radiation  is  restricted  to  a  range  between  -0.01  and 

48 
+0.24.   By  an  extension   of  the  formalism  for  intermediate  unobserved 

cascades  given  in  Appendix  1  it  is  possible  to  calculate  the  angular 
distribution  coefficients  for  the  9  -*■  1  radiation  which  result  from  the 
R  ->  19 ->■ 9  — ►  1  cascade .   It  is  found  that,  for  a  1+  resonance  level, 
the  minimum  possible  value  of  A„/An  occurs  for  an  assumed  spin  of  2  for 
level  19.   This  minimum  value  has  been  calculated  to  be  -0.175.   Thus 
a  1+  resonance  level  leads  to  a  minimum  possible  value  for  A_/An  of 
about  -0.093.   Comparing  this  with  the  observed  value  of  -0.24  +  .03  for 
the  2.34-Mev  gamma  rays  1+  can  be  ruled  out  as  a  possible  spin-parity 
of  the  resonance  level. 

Having  narrowed  the  possible  resonance  spin  values  to  2+  or  2-, 
the  contour  plots  may  be  examined  for  the  2+  ->»  1+  and  2--^l+  transitions, 
applicable  to  the  7.00-Mev  radiation  to  the  ground  state.   In  the  first 
case,  channel-spin  mixing  may  occur  while  in  the  second  case  proton- 
partial-wave  mixing  is  permitted.   For  the  2+-^l+  case  Fig.  16a  shows 
the  contour  plot  with  the  experimental  value  of  A  /A„  for  the  7.00-Mev 
gamma  ray  delineated  thereon  by  cross-hatching.   Fig.  16b  is  a  similar 
plot  for  A,/An.   The  darker  shading  in  Fig.  16b  indicates  the  superposi- 
tion of  cross-hatching  on  both  plots  and  defines  the  ranges  and  combi- 
nations of  mixing  parameters  which  can  give  rise  to  the  observed  coeffi- 
cient values.   For  £„  £ssO  (i.e.  Ml  radiation),  £,  is  constrained  to  be 
less  than  ^2 tan  50  .   This  must  not  be  taken  too  seriously,  however, 
since  the  A, /Aft  contours  are  widely  spaced  at  the  top  and  bottom  boundaries 
of  the  permitted  region,  so  that  a  very  slight  increase  in  the  error  for 
the  experimental  A,/A„  value  would  allow  S ,    to  range  up  to  90  .   The  small 
region  at  £„  ~tan(~67  )  and  £,  ^5  tan  49   also  fits  the  observed  corre- 
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lation.   The  plots  for  the  2 — +■!+   case,  shown  in  Figs.  17a  and  17b,  are 
quite  similar  to  the  2+->,l+  case  and  show  that  there  is  no  restriction 
on  f-p  proton-wave  mixing  for  &  £>S  0  (El  radiation).   The  two  small  per- 
mitted regions  at  £  p^tan  67  may  be  ruled  out  since  they  require  an 
M2/E1  mixture.   It  is  seen  that  an  unlikely  but  small  M2/E1  mixture  with 
£  ^iiO.026  or  greater  is  required  for  this  transition  for  the  case  of 
resonance  formation  by  pure  p  waves.   As  will  be  shown,  the  2-  resonance 
level  assignment  is  also  less  probable  than  2+  based  on  mixing  ratios 
which  are  required  by  experimentally  observed  coefficients  for  the  4.50- 
Mev  radiation. 

As  a  result  of  the  above  considerations,  a  spin-parity  of  2+  is  pro= 
posed  for  the  resonance  level  with  2-  as  a  less  likely  possibility. 

C.  R->13->0  cascade 

There  are  only  two  cascades  which  allow  cross  checking  of  mixing 
parameters  as  a  method  of  eliminating  possible  spin  sequences.   These 
are  the  cascade  through  level  19  and  that  through  level  13  both  of  which 
levels  have  undetermined  spins.   The  latter  cascade  will  be  considered 

first.   The  existence  of  an  A. /A~   term  which  has  been  observed  in  the  dis- 

40  23 

tribution  of  the  4.14-Mev  radiation  from  level  13  at  the  1686-kev  resonance 

is  evidence  that  level  13  must  have  a  spin  greater  than  or  equal  to  2.   Of 
the  likely  possible  spins  for  this  level,  2+s  3+,  and  4-,  only  2+  and  2- 
can  result  in  the  negative  value  of  A~/An  observed  for  the  4.14-Mev  gamma 
ray.   An  examination  of  the  contour  plots  for  the  four  sequences  2+-^2+->l+ 
for  the  R  ->13-^0  cascade  shows  compatibility  in  all  cases  with  the  ob- 
served coefficients  for  the  2.86-Mev  gamma  ray  and  the  4.14-Mev  gamma  ray. 
For  the  sequence  2+-^2--^l+  the  parameter  £   is  constrained  to  be  larger 
than  tan  70  but  the  other  possible  sequences  are  compatible  with  a  full 
range  of  values  for  either  £1  or  £.  as  the  case  may  be.   The  only  finding 
is  that,  for  all  cases,  the  R->13  transition  must  be  relatively  pure  di- 
pole  radiation.   The  observed  coefficient  values  for  the  weak  2.69-Mev  gamma 
ray,  which  also  is  believed  to  deexcite  level  13,  are  also  compatible  with 
all  four  possible  spin  sequences  with  reasonable  choices  for  o\  • 

D.  R->19->0  cascade 

Since  no  requirements  on  parameters  applicable  to  this  cascade  are 
imposed  by  transitions  already  considered  it  becomes  necessary  to  consider 
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the  R  ->19  ->0  sequence  with  no  prior  restrictions  on  the  deltas  involved, 
We  find,  furthermore,  that  no  information  is  forthcoming  from  the  R— >-19, 
2.50-Mev  gamma-ray  correlation  coefficient  since  this  radiation  is  indis- 
tinguishable from  the  12-^3,  2.48-Mev  gamma  ray. 

A  survey  of  maxima  and  minima  of  A9/An  coefficients  shows,  if  we 
allow  2-  as  a  possible  resonance  level  spin,  that  the  following  sequences 
for  the  R-^19-^0  transition  are  consistent  with  the  experimental  co- 
efficients for  the  4.50-Mev  radiation: 

2+  1+  1+  2+  1-  1+ 

2-  1+  1+  2-  1-  1+ 

2+  2+  1+  2-  2-  1+ 

2-  2+  1+  2-  3-  1+ 

2+  3+  1+  2+  3-  1+ 

It  has  been  found  in  the  analysis  that  a  better  fit  (see  Fig.  14) 

to  the  data  is  given  by  the  series  including  a  P,  term.   Then,  provided 

this  term  is  not  spurious,  which  is  not  to  be  excluded  as  a  possibility, 

the  list  can  be  reduced  by  eliminating  sequences  not  giving  rise  to 

A. /A_  terms. 
4  0 

The  following  possibilities  then  remain. 

2+  2+  1+ 

2=  2+  1+ 

2+  3+  1+ 

2+  3-  1+ 

2-  3-  1+ 

The  last  two  transitions  shown  require  E3/M2  radiation  so,  adhering  to 
the  hypothesized  El  and  E2/M1  radiative  transitions  as  practical  candi- 
dates, the  spin  of  the  4.50-Mev  level  is  limited  to  either  2+  or  3+.   Of 
the  three  remaining  possibilities,  the  2-  -^2+-Vl+  sequence  is  perhaps 
least  likely.   Figs.  19a  and  19b  show  that  this  sequence  for  R-^19-^0 
requires  an  f  to  p  wave  mixing  parameter  of  over  0.10  with  a  large  E2/M1 
mixing  ratio  of  nearly  12.0.   For  a  2+-^3+->l+  sequence,  Figs.  20a  and 
20b  show  that  an  E2/M1  ratio  of  over  1.7  is  required  for  the  4.50-Mev 
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radiation  in  order  to  fit  the  data.   Because  large  E2/M1  mixing  para- 
meters  are  not  commonly  seen  in  the  lighter  nuclei,  this  sequence,  though 
more  likely  than  2-  -*■  2+  -?»1+,  does  not  appear  as  probable  as  a  2+  •^•2+  -^-1+ 
sequence,  as  will  be  seen.   For  the  2+-^2+-^l+  sequence  Figs.  18a  and 
18b  show  that  for  an  assumed  channel-spin  mixing  ratio  of  zero,  an  E2/M1 
mixing  ratio  of  0.27  will  fit  the  observed  coefficients  for  the  4.50-Mev 
gamma  ray  with  a  considerable  latitude  permitted  for  Q  . .   Alternate 
choices  for  £,  will  also  permit  acceptable  choices  of  mixing  parameters 
for  this  sequence  though  in  general  &  and  Q-  will  be  required  to  be 
larger. 

Summarizing,  level  19  is  found  to  have  a  spin  of  2+  or  3+  under  the 
assumptions  which  have  been  made,  with  2+  being  more  probable. 

E.  Comparison  with  other  results 

Concurrently  with  the  present  work,  an  investigation  of  triple 

angular  correlations  at  the  1470-kev  resonance  has  been  undertaken  at 

40 
Wright-Patterson  Aeronautical  Laboratories.    The  angular  distribution 

of  the  7.00-Mev  gamma  ray  was  also  examined  in  this  experiment.   An  A_/A„ 

value  of  -0.38+.02  was  observed  for  this  gamma  ray  in  good  agreement  with 

the  value  shown  in  Table  7.   Preliminary  results  of  this  experiment  have 

indicated  restriction  of  the  resonance  level  spin  to  a  value  of  1  or  2 

and  show  that  a  zero  spin  for  the  4.50-Mev  level  is  not  compatible  with 

such  a  restriction.    The  preliminary  analysis  of  the  triple-correlation 

experiment  is  seen  to  be  in  agreement  with  the  present  work.   It  is  to  be 

hoped  that  the  complete  analysis  of  this  work  which  is  now  in  progress, 

will  lead  to  a  firmer  indication  of  the  spin-parity  of  the  4.50-Mev  level 

since  the  2+  (3+)  value  proposed  herein  depends  upon  the  somewhat  debatable 

A,/An  term  observed  and  upon  the  hypothesis  that  only  El  or  E2/M1  radiation 

is  present. 

F.  T  =  1  character  of  level  19 

The  prediction  of  0+  for  the  spin-parity  of  level  19  is  based  on  the 

30 
premise  that  this  level  is  the  T  =  1  analog  of  the  3.80-Mev  level  of  Si 

42 
in  the  isotopic  spin  multiplet  formalism.    A  spin  of  0  has  been  suggested 

30 
for  the  3.80-Mev  level  in  Si   based  on  observations  of  cascade  gamma  rays 

following  the  Si   (p,p")Si   reaction.    A  spin  of  0+  results  also  from 
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calculations  made  by  Thankappen  and  Pandya  of  collective  vibrational  states 

in  a  weakly-coupled  collective  model  of  the  nucleus.    The  spin-parity  of 

30 
level  19  in  P   is  thus  expected  both  on  experimental  and  theoretical  grounds. 

The  finding  of  non  0+  spin-parity  for  this  level  then  appears  to  be  good  evi- 
dence that  either  the  level  is  complex  (e.g.  a  doublet)  or  that  it  is  not  the 

30 
isotopic  spin  multiplet  analog  of  the  4.80-Mev  state  in  Si   as  has  been  ac- 

30 
cepted.   Some  evidence  has  been  reported  that  the  3,80-Mev  level  in  Si   is 

a  doublet,  based  upon  observations  of  protons  from  the  reactions 

0.28,   .„.30    .  „.30,    ,N„.30   59,44  _   .      ...   t,    ,      . 
Si   (t,p)Si   and  Si   (P,P  )Si   .         It  is  possible,  therefore,  that 

30 
the  analogous  4.50-Mev  level  in  P   is  a  doublet,  although  no  other  evi- 
dence so  suggests.   The  other  possible  explanation  is  that  a  nearby  level 

30 
in  P   ,  such  as  that  at  4.42  Mev,  may  be  the  T  =  1  level  rather  that  the 

4.50  Mev  level.   This  does  not  appear  unlikely,  since  the  T  =  1  character 

attributed  to  the  4.50-Mev  level  is  a  probable  rather  than  a  firm  assignment 

32      30 
which  was  deduced  from  the  investigation  of  the  S   (d,a)P   reaction  by  Endt 

18  30 

and  Paris.     It  is  based  upon  a  correction  of  the  3.80-Mev  level  in  Si 

for  the  Coulomb  energy  and  neutron-proton  mass  difference  which  yields  a  pre- 

30 
dieted  energy  of  4.47  Mev  for  the  excitation  energy  in  the  P   nuclide.   The 

assignment  is  strengthened  by  the  fact  that  the  alpha-particle  group  from 

32      30  30 

S   (d,CC)P   corresponding  to  a  residual  excitation  in  P   of  4.501+  .010  Mev 

has  only  about  half  the  intensity  of  neighboring  groups.   This  has  been  inter= 

preted  as  being  due  to  operation  of  the  selection  rule  conserving  isotopic 

spin  in  this  reaction. 

30 
On  the  basis  of  the  energy  predicted  for  the  T  =  1  level  in  P   the 

4.42-Mev  level  is  seen  to  be  nearly  as  likely  a  candidate  as  the  4.50-Mev 

level,  particularly  since  the  energy  prediction  is  not  expected  to  be  exact. 

Similarly  in  the  T  conservation  argument,  the  observed  lower  intensity  of 

the  alpha  group  for  the  4.50-Mev  level  can  hardly  be  considered  conclusive 

evidence  of  the  T  =  1  character  of  the  level,  since  some  T  =  0  levels  are 

associated  with  alpha-particle  groups  of  even  lower  intensity  (e.g.  the 

30 
3.84-Mev  level).   Unfortunately,  the  4.42-Mev  level  in  P   is  not  known  to 

be  excited  in  radiative  capture,  so  that  spin-parity  of  this  level  will  pro- 
bable have  to  be  determined  from  a  different  reaction. 
5.   Summary  of  Results 

Excitation  curves  for  the  Si   (p, 5  )P   reaction  have  been  obtained  over 
the  proton  energy  region  from  1420  through  2160  kev.   Four  previously  unreported 
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resonances  have  been  detected  and  confirmed  at  proton  energies  of  1975, 
2033,  2075,  and  2117  kev  with  experimental  uncertainties  of  +4  kev. 
Gamma-ray  spectra  at  these  resonances  show  that  inelastic  scattering  pre- 
dominates  over  radiative  capture,  except  at  2117  kev  where  no  inelastic 
scattering  is  evident.   Some  information  concerning  the  radiative  decay 
schemes  at  the  1975,  2075,  and  2117  kev  resonances  was  obtained  by  spectral 
analysis  with  the  assistance  of  sum-coincidence  spectra. 

The  gamma-ray  spectrum  at  the  1470-kev  resonance  has  been  analyzed 

23 
and  found  to  be  in  agreement  with  the  decay  scheme  proposed  by  Moore. 

The  angular  distributions  of  gamma  rays  observed  at  this  resonance  have 
been  measured  in  a  series  of  five  experimental  observations,  each  con- 
sisting of  gamma-ray  spectra  obtained  at  five  angles  relative  to  the  pro- 

ton  beam.   Analysis  of  this  data,  under  the  assumption  that  El  and  E2/M1 

30 
radiation  deexcite  P   energy  levels  of  several  Mev,  has  shown  that  the 

1470-kev  resonance  level  has  a  spin-parity  of  2+  or,  less  probably,  2-. 

The  spin-parity  of  level  13  at  4.14  Mev  has  also  been  found  to  be  limited 

to  2+  or  2-.   Level  19  at  4.50  Mev  was  found  not  to  have  a  spin-parity  of 

0+  which  had  been  predicted.    Under  the  above-mentioned  assumption,  an 

apparent  A./A-.  term  observed  in  the  angular  distribution  of  this  radiation 

leads  to  the  conclusion  that  the  spin-parity  for  level  19  is  2+  or,  less 

probably,  3+. 
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VI  APPENDICES 
Appendix  1  -  Theoretical  Description  of  Double  Angular  Correlations 

A .  General 

There  has  been  abundant  literature  coverage  of  the  theory  of  angular 
correlations  of  successive  nuclear  radiations  since  Hamilton  first  pre- 
dicted the  effect  and  developed  theoretical  expressions  for  successive 
emission  of  quanta,  ^  A  comprehensive  review  of  the  subject  has  been 
published  by  Biedenharn  and  Rose  ^  and  recent  articles  useful  in  exper- 
imental applications  have  been  published  by  Gove  ^'  and  by  Biedenharn*  ^° 
The  following  notation  is  essentially  that  of  Biedenharn.  In  all  cases  a 
primed  symbol  represents  a  larger  value  than  the  same  symbol  without  a 
prime. 

In  the  case  of  a  (p,^*)  reaction  in  which  the  direction  of  a  primary 
gamma  ray  deexciting  a  resonance  level  is  observed  relative  to  the  direc- 
tion of  an  incident  proton  beam^  the  angular  correlation  is  referred  to 
as  a  double  correlation.  The  function  describing  the  relative  probability 
of  the  gamma  decay  as  a  function  of  the  single  between  the  beam  direction 
and  the  emitted  gamma  ray  can  be  written  as  an  expression  in  the  even 
Legendre  polynomials % 

Z^max 
W(e)  s  2   Ajl)A  (2)P^(cose)       ,   even       (1) 

u  r  0 

where   A„(i)  Z  ajL^Fjh^^J)  +  26^^^)?^!^^) 

+  6ha*<LiLX<LlL^iJ> 
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1 
2  . 


and  Fj.CLL'jiJ)  I  (-l)^""1"1  [(2L+1)  (2L+l)(2J+l)J  '< 

(LL81  -1|V0)W(JJLL!;Z/Ji) 

The  index  h  is  related  to  the  index  i  by: 

h  r  i2  -  i  +  1 

(LL  1  -l|Z/0)  and  W( JJLL  J  V  ^)     are  respectively  the 

Clebsch-Gordon  coefficients  and  the  Racah  coefficients  which  are  commonly 
used  to  describe  systems  involving  coupling  of  angular  momenta.  The  F 
functions,  for  which  an  extensive  tabulation  is  available,    are  defined 
in  terms  of  a  Clebsch-Gordon  and  a  Racah  coefficient  by  Biedenharn  and 
Rose.  ^   Arguments  of  the  coefficients  are  angular  momenta  of  the  se- 
quence: 


(2) 
2} 

j1      =  channel  spin  s  proton  spin©  target  spin 

L^  ,1*!     =  interfering  proton  orbital  angular  momenta 

j     =  spin  of  resonance  level 

L  L1     =  interfering  gamma  ray  mult ipolari ties 
2*  2 

jp     =  spin  of  final  state  (not  necessarily  a  ground  state) 
L      -  mixing  parameter  representing  the  ratio  of  the  reduced 
matrix  element  for  the  transition  involving  LJ  proton 
waves  to  that  for  L^  waves 

£      =  mixing  parameter  representing  the  ratio  of  the  amplitude 

i 

of  gamma  rays  of  raultipolarity  L2  to  those  of  multipo- 

larity  L2 
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a(LL')  =  co4L-|L.)  -  2[l(L+1)l'(lViJJ* (3) 

is  a  particle  parameter  in  which  c,  and  &Ttare  the  phases  at  the  nuclear 
surface  of  the  particle  partial  wave  functions  of  multipolarity  L  and  L  . 
The  coherent  interference  between  partial  waves  is  expressed  by  the 
cos(£t  "^l')  term.  The  factor  «  includes  the  coulomb  phase  shift  due  to 
electrostatic  potential  for  charged  particles  and  the  phase  shift  occur- 
ring in  "hard  sphere"  scattering.  The  si-CL^o)  terms  for  photons  as 

opposed  to  particles  are  equal  to  unity.  In  the  case  where  mixing  of 

i 

channel  spins  j^  and  j^   occur,  W(6)  may  be  expressed  as  an  incoherent 

mixture : 

w(e)  =  wj;L(e)+  $2  w^(9)  (4) 

r2 

where  q,  is  a  mixing  parameter  expressing  the  ratio  of  probability  of 

channel  spin  jn,  entering  the  reaction  to  that  of  channel  spin  Ji,(Jt  j-j). 

The  maximum  order  of  Legendre  polynomials  appearing  in  (l),  ZJ         ,   is 

max 

determined  by  conditions  imposed  on  the  W  (abed;  ef)  coefficients  by  con- 
servation of  angular  momentum  which  may  be  expressed  by  the  so-called 
triad  relationships.  These  require  each  of  the  triads  of  values  (acf), 
(bdf),  (abe),  and  (cde)  to  form  a  possible  triangle.  Restriction  to  the 
case  of  intermediate  states  of  well  defined  parity  (i.e.  no  interf erring 
intermediate  states  of  opposite  parities)  leads  to  the  further  simplifi- 
cation that  U   may  take  on  only  even  values. 

In  the  case  of  the  reaction  (p,x^),  in  which  x  represents  an  unob- 
served intermediate  gamma  ray,  the  correlation  is  expressed  as  an  incoher- 
ent mixture  of  multipoles  occurring  in  the  unobserved  transition,  i.e.  in 
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the  sequence: 


(5) 


r  2 

W(e)  ctn  b3  expressed?     W(0)     Z    \ ,(9)  +•  O  2Wti  (6)  (6) 

where: 

jQ  a  channel  spin 

Iq,Lq  -  interfering  proton  orbital  angular  momenta 

J,  a  spin  of  resonance  level 

LpL^  a  interfering  gamma  ray  multipolarities  of  unobserved 

radiation 

J  a  spin  of  intermediate  level 

Lp,L2  -  interfering  gamma  ray  multipolarities  of  observed  radiation 

j   =  spin  of  final  state 
2 

<^P  -  ratio  of  intensity  of  the  unobserved  gamma  ray  of  higher 

multipolarity  to  the  intensity  of  the  lower  multipolarity 

gamma  ray. 

The  correlation  function  now  contains  an  additional  part  expressing 

the  connection  between  the  two  intermediate  levels  by  the  unobserved 

radiation  as  follows: 

Z'max 
W(6)  ^  A1>(0)Cs>(l)Aj2)PE>(cos9)  (7) 


where  the  A  (i)  are  the  same  as  in  (1)  and 

Cv(l)     Z     (  -=1)V[(2J11-1)(2J2-I-1)]^(J1     LjJ2*  J1J2> 

It  is  seen  that  for  the  general  double  correlation  case  four  mixing 
parameters  may  occur  (barring  the  further  complication  of  mixed  inter- 
mediate states) ♦ 
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In  the  present  investigation  the  Si^9  target  nucleus  has  a  ground- 
state  spin  of  l/2 +,  with  the  simplifying  result  that  channel-spin 
mixing  and  proton-partial-wave  mixing  may  not  occur  simultanously0  The 
number  of  6 's  is  thus  limited  to  a  maximum  of  three.  Explicitly ^  for 
resonance  levels  of  1+-  ,2-,  3\»  U->   etc  orbital  angular  momentum  mixing 
of  two  partial  waves  may  occur  but  channel-spin  mixing  is  prohibited,, 
For  the  resonance  level  spins  of  opposite  parities  from  those  listed 
only  channel-spin  mixing  is  permitted. 

When  three  mixing  parameters  occur  in  the  general  case  of  inter- 
mediate unobserved  radiation  (p,xtf )  with  a  spin  sequence 

or    JnlNll   kl  2|J2  (8) 

\iol 

the  normalized  correlation  function  can  be  expressed  as 

w(e)  =  i/AqZ  A^(cose)  (9) 

where  each  A  /A  is  a  product  of  the  forms 


AxA>  =  TT^(^)  (105 

withy    =     1,2,3  or  2,3,4 


and  g       z 

(It  is  to  be  noted  that  the  coefficient  A  v  appearing  in  (9)  is  a  conven= 
tional  representation  of  the  more  complete  forms  shown  in  (1)  and  (?). 
A  distinction  must  also  be  made,  of  course,  between  the  correlation  func= 
tion  W(0)  and  the  Racah  coefficient  W(abcd;  ef). 
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The  coefficients  a,  b,   and  c  can  te  expressed  as  follows: 


/  a 


protcn 

wave 

mixing 


*1 


2Lo(L0+  1)   -2>^+l) 


■U 


^1 


s     4[Lo(LQ  4  1)1^  ±  1)1 »  FjLqL^) 


Lo(lo  +1)+  L£(l£  +  1)  -  *>(*>  + 


2L6(Lq  j;  1) 


2Lo(Lq-»-  1)  -  M^+D 


F»»(IioJoJi^ 


:u: 


unobs0 
gamma      ^ 
mixing 


\/2 


°*2 


c«/2 


(-if  [(2J!*  1)(2J2+  I)}*  W(«W23  JjJg) 


(-1)X/[(2J1+  1)(2J2  +  1)]*  WCJj^LJJ^Jj 


gamma 
mixing 


a    o 


<   b~3 


c*3 


F^(L2^2J2) 
2FjL2L2j2J2) 

F„a2J2J2) 


f  a,4    -    2^0. 


+  i: 


2Lo(L0  +  1)  -  *(*>  + 


chan, 

spin       \    ^4    " 

mixing 


2Lo(Ln  +  1) 


F*,(%Vl^ 


*"*  2L0(Lo+  1)  -v(vU)  °0J 


t2    =    1\    "     *"?4    -    : 
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proton 
wave 

mixing 


For  the  (p,  8"  )  case  the  only  changes  required  are:  &iid~  ~9   \2  S  ®» 

c^2=  °  and  Ji  =  J2  =  J° 

The  angular  correlation  function,  as  an  alternative s   may  be  ex- 
pressed in  terms  of  the  coefficients  Z  and  Z^  9   instead  of  the  F  coef- 
ficients* 

b^  =  2(-l)J°Z(l0J1I^J1,  Jo")  (12)' 


,30 


»  .  -  « 


c*i  =  (-1)    «WiVi»  v> 


c  a 


unobs. 

gamma 

mixing 


u2    =  (-1)°  X  W( J^gJ^g;  Iq.  2/) 


D*2 
^  c«,2 


r  0 


(-i)°Li  wC^J^Jgj  l[i>) 


gamma 
mixing 


a 


v   c^3     s 


(-l)^2  Z1(L2J2L2J2,   j2i/) 
2(-l)J2  Z^J^   j2v) 
(-1)^2  ^(L^J^Jgi  j2v) 


/    a 


Jo 


chan. 

spin 

mixing 


iw4     =     (-1)  U  ZCLqJ-lLqJ-l,   j02^) 


^4 


=     (-l)Ji  ZCLqJ^o^i   JqV) 
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For  (p,  2f)  reactions  J.  -  Jp  ■  J,  a^  2  *  (2J+1)~  ,  b^g  s  0S   c^  s  ° 

*  (  The  two  formalisms  (11)  and  (12)  are  to  some  extent  incompatible . 

.   50  51 
Not  only  may  the  coherent  mixing  parameter  phases  differ  by  II  '  *' 

but  also  for  L =  0  the  parameter  S,  doesn't  represent  the  same  matrix 

element  ratio  in  both  notations.  The  notation  above,  (12),  has  been  used 

for  most  calculations  in  the  present  work.  These  expressions  are  in 

agreement  with  Sharp  -*  except  that  the  sign  of  bv  o  ,  for  the  (p,xfl) 

case  only,  has  been  changed  to  agree  with  the  part  of  Huby's  correction 

applying  to  the  observed  ganima  ray.'*  *'   For  the  "unnatural"  radiation 

mixtures  (ML/E2,  23/^2*  etc.)  the  sign  of  fcuj  should  be  reversed. 

This  sign  convention  penaits  the  use  of  Moore's  tabulated  values  2-* 

where  the  sane  convention  is  used.) 

The  Z  and  Z,  coefficients  and  Racah  and  Clebsch-Gordon  coefficients 

52  °^ 

have  been  tabulated  by  Sharp  et  al.     Using  this  tabulation  Moore  *'-' 

has  prepared  a  convenient  table  of  the  a,  b,  and  c  coefficients  applica- 

29 
ble  to  radiative  proton  capture  by  Si  . 

In  order  to  fit  experimental  data  on  angular  correlation  coeffi- 
cients to  theory,  it  is  desirable  to  first  eliminate  as  m?ny  spin  se= 

quences  as  possible  on  the  basis  of  any  evidence  which  may  be  available. 
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In  cases  where  only  two  mixing  parameters  enter,  contour  plots 
ing  contours  of  constant  theoretical  values  of  the  coefficients  A2/Aq 
can  be  plotted  against  the  inverse  tangents  of  the  two  mixing  parameters 
for  all  possible  spin  sequences »  If  three  mixing  parameters  enter $   one 
of  them  may  be  held  constant  for  the  purpose  of  plotting  the  correlation 
coefficient  contours  against  the  inverse  tangents  of  the  other  two  param- 
eters It  is  seen  from  (10)  that  the  contour  values  for  different 
values  of  the  fixed  parameter  differ  only  by  a  multiplicative  constant. 
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When  experimental  values  of  correlation  coefficients  are  delineated  on 
the  contour  plot,  some  spin  sequences  can  be  eliminated  and  the  possible 
combinations  of  mixing  parameters  are  restricted  for  other  sequences. 
The  observed  values  of  several  correlation  coefficients  in  one  decay 
scheme  may  then  fix  or  highly  restrict  values  of  the  spin-parities  of 
resonance  and  intermediate  levels,  and  determine  the  values  of  the  mix- 
ing parameters. 
B.  Application  to  1470-kev  resonance 

For  the  resonance  occurring  at  a  laboratory  proton  energy  of  1470 
kev,  the  phase  factor  cos'£-  -|t()  has  been  calculated  using  the 
expression^  =^l+,^l,  wnere  <J.  represents  the  phase  shift  due  to  the 
point  coulomb  potential  andp^  is  the  "hard  sphere"  phase  shift. 

tfT  is  defined  k2   by  exp(2i<ST  )  s  P(L  +  1+  1*[  ) 

L  L       ra  +  i-m) 

22*  e2 

where  ns   ■  which  leads  tos 

The  value  for  <£   defined  by 

exp(2ifL)  -  GL(^kR)^iFL(TkR)      ^ 
GL(^,kR)  H-iFL(^,kR) 
tfhere  F  and  G  are  the  regular  and  irregular  solutions  to  the  radial  wave 

equation  for  the  coulomb  potential.  The  value  for  (L     can  be  found  from 

53 
the  graphs  of  Sharp  et  al.     (these  graphs  actually  list-d&  )  and  em- 

1/3 
ploying  R  »  1.5A  '  J£„   as  an  estimate  of  the  nuclear  radius.  These 

graphs  have  also  been  used  to  calculate  the  relative  pentration  factors 

vlal-2  =  %*:%?  (i6) 

which  are  useful  as  a  guide  to  the  expected  mixing  parameter  of  inter= 
fering  partial  proton  waves.    Results  of  these  calculations  appear  in 
Tables  8  and  9. 
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TABLE  8 
Relative  Penetration  Factors  for  1470-Kev  Protons 


d/s 


f/P 


g/d 


3.1  x  10 


-2 


5o0  x  10 


1.3  x  10 


TABLE  9 
Coulomb  Phase  Factors  for  1470-Kev  Protons 


*,-*. 


104 


3   Si 


73* 


4   T2 


56 
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Appendix  2  -  Computer  Programs  and  Analysis 

Ao  General 

In  order  to  analyze  experimental  data  in  the  form  of  gamma-ray 
spectra  recorded  by  the  128-channel  analyzer,  two  computer  programs 
were  developed.   The  first  and  more  extensive  program  entitled  program 
SPECFIT  computes  experimental  correlation  coef f icientss  while  program 
CONTOUR  is  designed  to  utilize  a  California  Computer  Products  Model 
565  graph  plotter  as  an  output  medium  to  automatically  plot  contours 
of  constant  angular  correlation  coefficient  values  against  the  inverse 
tangents  of  two  mixing  parameters .   Both  programs  are  written  in  fortran 
language  and  were  run  on  a  CDC-1604  digital  computer. 

B.   Program  SPECFIT 

The  output  pulses  of  a  photomultiplier  tube  viewing  a  scintillation 
crystal  upon  which  monoenergetic  photons  impinge  do  not  all  have  the  same 
voltage.   Instead,  the  three  fundamental  processes  of  photoelectric  ab- 
sorption, Compton  scattering  and  pair  production,  which  occur  in  the 
crystal  with  probabilities  dependent  upon  gamma-ray  energy,  lead  to  a 
response  pulse-spectrum  which  is  characteristic  of  each  particular  gamma- 
ray  energy.   Experimental  geometry,  crystal  size,  shielding,  various  line- 
broadening  effects,  and  characteristics  of  individual  crystal-photomulti- 
plier  tube  combinations  influence  the  response  also,  so  that  response 
functions  are  in  general  unique  for  each  spectrometer  in  a  given  experi- 
mental arrangement. 

An  experimentally  observed  spectrum  is,  of  course,  usually  composed 
of  a  number  of. gamma  rays,  the  spectrum  being  the  sum  of  the  constituent 
radiations  each  multiplied  by  its  intensity.  The  problem  of  deciphering 
such  an  observed  spectrum  to  obtain  the  energies  and  intensities  of  com- 
ponent gamma  rays  is  of  fundamental  importance  in  gamma  ray  spectroscopy 
and  has  not  yet  been  satisfactorily  solved  in  complete  generality. 

Using  notation,  which  is  appropriate  for  the  discrete  pulse  height 
intervals  into  which  multi-channel  analyzers  sort  pulses,  and  also  ap- 
propriate for  the  discrete  number  of  gamma  ray  energies  composing  re- 
sonance deexcitation  spectra,  one  can  write; 

m 
S_.  =  " 


J  =  l 
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where  S.  =  observed  number  of  counts  in  channel  i 

1 


R. .  =  probability  that  a  photon  of  energy  E. 
originating  at  a  particular  source 


will  produce,  a  pulse  falling  in  pulse- 
height  channel  i  (usually  referred  to 
as  the  response  matrix) 


I.  =  the  number  of  photons  of  energy  E. 
**    originating  at  the  source 

m  =  the  number  of  gamma  ray  energies  which 

occur  in  the  spectrum  under  consideration 


A  rather  fundamental  and  obvious  approach  to  this  problem  is  to 
assume  no  knowledge  of  gamma-ray  energies  is  available,  divide  the  gamma 
ray  energy  range  into  a  number  of  discrete  intervals  equal  to  the  number 


of  channels,  then  solve  for: 


i  '   ^-r1  4-»   ik   Kj  j    £+       ik   k 


In  fact,  the  first  attempt  to  unfold  experimental  spectra  in  the 
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present  work  utilized  a  modification  of  a  program   which  employed  an 

iterative  technique  developed  by  Schofield  to  solve  the  above  equation. 
This  attempt  failed  for  two  reasons  which  are  rather  instructive.  First 3 
the  analytical  expression  used  to  approximate  the  response  matrix  did  not 
adequately  represent  the  true  spectrometer  response,  and  second,  the 
technique  proved  inadequate,  with  experimental  counting  statistics  attain- 
able to  handle  the  analysis  of  approximately  100  channels. 

In  order  to  avoid  the  first  problem  resort  was  had  to  the  technique 
of  interpolating  between  experimentally-determined  monoenergetic  cali- 
bration curves  in  order  to  obtain  the  response  matrix.   In  order  to  over- 
come  the  second  objection  it  was  decided  to  use  information  on  the  gamma 
ray  energies  available  from  preliminary  "stripping"  of  spectra  and  from 
other  sources.   Having  thus  reduced  the  number  of  gamma  rays  to  be  con- 
sidered, a  linear  combination  of  the  selected  response  matrix  elements 
was  then  fitted  to  the  observed  spectrum  by  a  least  squares  process. 
This  procedure  was  followed  by  another  least  squares  analysis  to  obtain 
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correlation  function  coefficients  from  the  gamma  ray  intensities  cal- 
culated for  a  series  of  runs  at  different  experimental  angles.   The 
steps  in  the  computation  procedure  in  program  SPECFIT  as  finally  de- 
veloped are  as  follows: 
L   Load  data.   (See  section  D  of  this  Appendix) 

2.  Subtract  any  constant  background  from  calibration  spectra.   This 
background  represents  the  Compton  "tails"  of  higher  energy  gamma  rays 
contaminating  certain  of  the  calibration  gamma  rays. 

3.  Normalize  each  calibration  spectrum  to  a  photopeak  height  of  one 
unit . 

4.  Include  a  straight  line  in  the  response  matrix  if  this  option  is 
selected.   This  step  is  essentially  a  refinement  of  step  2. 

5.  Correct  the  input  and  background  spectra  for  dead  time  of  the  pulse 
height  analyzer.   For  this  purpose  the  dead  time  given  by  the  manufacturer 
is  38  +  0.5i  microseconds  where  i  is  the  pulse  height  channel  number. 
This  leads  to  an  expression  for  corrected  counts; 

N. 

N'.  =  -         -L 


j  127 

1  -  1/T  2T  N-(38  +  0.5i)»10"6 
1^2   L 


where   N.  =  the  number  of  counts  observed  in  channel  i 

l 

T  =  counting  time  in  seconds 

6.  Calculate  energies  in  Mev  corresponding  to  each  channel  number  of 
the  input  spectrum. 

7.  Adjust  energy  scale  and  number  of  channels  of  background  spectrum 
to  match  the  input  spectrum.   Multiply  by  ratio  of  run  times  to  get  the 
proper  background  to  be  subtracted  from  the  input  spectrum  and  perform 
the  subtraction. 

8.  Using  the  monoenergetic  calibration  spectra,  calculate  a  matrix  of 
response  curves  for  the  gamma  rays  which  the  spectrum  is  to  be  decomposed 
into.   This  calculation  is  performed  in  a  subroutine  entitled  RESPONS  as 
follows : 

a.   Find  the  two  calibration  curves  bracketing  a  gamma  ray  energy 
desired  as  a  response  matrix  element. 


75 


bo   Shift  the  lower-energy  calibration  curve  found  in  a.  so  that 
its  photopeak  falls  in  the  same  channel  as  the  higher  energy 
photopeak.   Fill  in  bottom  of  shifted  curve  as  determined  by 
the  slope  between  the  bottom  two  unshifted  channels.   Match 
energy  calibration  of  the  lower-energy  curve  to  the  upper. 
(For  this  step,  and  elsewhere  where  interpolation  is  necessarys 
cubic  interpolation  is  performed  using  the  four  nearest  points 
to  the  desired  value,) 

Co   Interpolate  linearly  between  each  pair  of  points  of  the 

bracketing  calibration  curves  to  obtain  a  response  curve  for 
the  desired  gamma-ray  energy.   (A  "curve"  in  this  description 
is  meant  to  imply  the  vector  consisting  of  a  1  X  n  array  of 
values,  n  being  the  number  of  channels  used  in  the  description). 

d.  Using  the  input  spectrum  channel  energies  calculated  in  6.s 
interpolate  to  get  a  response  curve  with  channels  and  channel 
energies  corresponding  to  those  of  the  input  data. 

e.  Normalize  the  response  curve  for  the  given  gamma-ray  energy  so 
that  the  ordinate  corresponding  to  each  channel  represents  the 
probability  that  one  photon  leaving  the  source  will  produce  a 
count  in  that  channel.   For  this  calculation  the  experimentally- 
determined  intrinsic  peak  efficiency  data  of  Lazer   h«ve  been 
used  to  derive  the  following  approximate  expression  for  intrinsic 
peak  efficiency  as  a  function  of  source  distance: 

€  =  AE"X 
where   A  =  . 165  +  .0067h 

x  =»  .9699  -  .00918h 

h  =  source  to  crystal  distance  in  centimeters 

E  =  gamma  ray  energy  in  Mev 

Assuming  an  accurately  gaussian  photopeak,  and  that  the  crystal 
resolution  is  approximately  inversely  proportional  to  the  square 
root  of  gamma  ray  energy,  (an  assumption  which  was  experimentally 
verified)  the  ratio  of  photopeak  height  to  area  can  be  expressed: 


"   «662*E 


^r_  WR      f    „bbZjtE  866wr    j_ 
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where  R  is  the  usual  resolution  definition  of  peak  width  at 

137 
one  half  maximum  amplitude  divided  by  energy  for  the  Cs 

gamma  ray,  and  W  is  the  spectrometer  calibration  in  pulse 
height  channels  per  Mev.   When  OC   is  multiplied  by  the  appro- 
priate solid  angle  fraction  and  attenuation  coefficient,,  the 
value  obtained  represents  the  probability  of  recording  one 
count  in  the  photopeak  channel  when  one  photon  originates  at 
the  source. 
f .   Return  to  main  routine. 
9o   Fit  the  experimental  spectrum  with  a  linear  combination  of  the  selected 
monoenergetic  response  curves,  the  coefficients  representing  absolute  in- 
tensities because  of  the  method  of  normalizing  the  response  matrix.   The 
weighted  sum  of  the  squares  of  the  differnces  between  observed  channel 
counts  and  fitted  counts  is  minimized  in  order  to  arrive  at  this  fit. 
The  normal  equations  resulting  from  the  minimization  are: 

2S.RMw.  =  V  R..R..I.W.  =  0    i<  j 
j  ij 


where  S.  =  Observed  counts  in  channel  i 
J 

R   -  Response  matrix  element  giving 
J    probability  of  a  count  in  channel 
j  due  to  one  photon  of  energy  E, 
originating  at  the  source 

w.  =  Weight  assigned  to  the  squared 
J   error  in  channel  j 

I.  =  Number  of  photons  of  energy  E. 
which  leave  the  source  during 
the  counting  period 

j  summations  are  over  the  channels  for 
which  a  fit  is  desired 


i  summation  is  over  the  number  of  gamma 
ray  energies  E.  used  to  fit  the  ob- 
served spectrum 


For  calculation  the  I.  are  expressed;   I.  =  J>      A,  .  B, 

x  l  *-f       ki  k 

k 

with  A.  .  =  J.     R..RMw.      and       B,  -  y°[    S.R..W. 
ki    ^   Ji  jk  J  k    4^   j  jk  j 
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A  standard  subroutine  (GAUSS3)  is  used  to  invert  A,   Mention  should  be 

made  of  the  weights  w.  which  have  been  used.   Normally ,  for  a  Poisson 

distribution  the  weight  is  taken  as  l/£   where  the  number  of  observed 

2 
counts  in  a  channel  are  used  as  an  estimate  of  the  variance  &  . 

Ferguson,    who  has  pointed  out  that  this  procedure  overweights  statist!- 

cally  low  counts,  identifies  the  Poisson  distribution  variance,  which 

rightly  is  the  mean  number  of  counts,  with  the  fitted  number  of  counts 

for  each  channel.   This  procedure  is  adopted  here  with  a  modification  to 

account  for  statistical  error  arising  from  the  subtraction  of  background 

counts.   We  have  S!  =  S.  =  b.,  where  S.  are  observed  channel  counts  cor- 

111         1 

rected  for  dead  time  and  b.  are  background  counts  to  be  subtracted.   Then,, 


<£2(s!)  =<£2(s.)  +  £2(b.) 


«*  S.  +  (T/T,  )b.  =  S!  +  (1  +  T/T,  )b. 
l    v    b   i     1   ^       W    l 

where  T  =  counting  time  for  data  and  T,  =  counting  time  for  background 
run.   Identifying  S1.  with  the  computed  number  of  counts,  S9.  '  ,  then  gives 

S! '  +  (1  +  T/T,)b.  for  the  channel  variances.   To  avoid  overweighting  of 

1  ^2 

very  small  values  of  fitted  counts  an  arbitrary  minimum  value  of  S   is 

^2 
chosen  to  be  two  per  cent  of  the  average  value  of£.  over  the  fitted 

channels.   Since  no  fitted  count  values  are  initially  available,  an  ini- 
tial weight  of  one  is  assigned  to  each  channel  and  the  solution  is  then 

iterated  to  convergence,  stopping  when 

2 

=  6 


Ft* 


is  smaller  than  10 


(This  series  of  iterations,  usually  five  in  the  present  work,  accounts 
for  most  of  the  program  running  time  and  can  be  abbreviated  for  prelimi- 
nary analyses;  see  reference  57  in  this  regard.)  CI. /I.  is  the  fractional 
change  in  the  computed  intensity  of  the  i  "  gamma  ray  between  successive 
iterations. 

10.  The  computed  intensities  are  examined  for  algebraic  sign.   Since 
only  positive  intensities  have  a  physical  basis  any  negative  intensities 
are  discarded  and  the  problem  resolved  to  this  point. 

11.  For  each  gamma  ray  calculate  the  fitted  total  number  of  counts  in  the 
channels  above  the  channel  representing  40  per  cent  of  the  photopeak  energy, 
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(This  is  an  arbitrary  calculation  used  only  for  checking  results  obtained 
in  manual  unfolding) .   The  fractional  error  assigned  to  these  counts  is 
set  equal  to  the  fractional  error  in  the  corresponding  intensity .   The 

statistical  variance  estimators  of  the  intensities  have  been  taken  to  be 

41 
the  diagonal  elements  of  the  inverse  of  the  matrix  A  in  9.     In  addition 

to  statistical  intensity  errors  there  is  an  error  due  to  gain  and  discrimi- 
nator drifts  in  the  monitor  counts  channel  which  must  be  accounted  for  as 
well  as  error  due  to  eccentricity  of  the  beam  spot  on  the  target.   These 
fractional  errors  have  been  estimated  to  be  3.3  per  cent  and  1.2  per  cent 
respectively  in  the  present  experiment.   The  three  errors  are  combined 
according  to; 

3 

Two  other  possible  sources  of  error  should  be  noted  here.   First,  the 
errors  in  the  experimentally  determined  response  matrix,  which  can  con- 
siderably complicate  the  error  calculations  for  the  intensities,  are  kept 
to  a  minimum  by  obtaining  good  counting  statistics  in  the  calibration 
spectra.   The  other  source  of  error  is  gain  drift  which  tends  to  smear 
spectral  peaks.   Fortunately,  the  main  effect  of  such  drift  seems  to  be 
an  increase  in  the  mean  squared  errors  in  the  curve  fitting  process  and 
intensity  determinations  are  apparently  not  too  sensitive  to  the  effect. 

12.  The  intensities  and  their  errors  calculated  in  11  are  read  out  and 
stored  for  use  in  the  subsequent  determination  of  correlation  coefficients 
and  mean  values. 

13.  The  absolute  values  of  differences  between  observed  and  fitted  counts 
are  summed  over  fitted  channels  and  divided  by  total  input  spectrum  counts 
over  the  same  range  to  give  a  measure  of  the  fit  which  has  been  achieved 
in  the  analysis.   This  measure  is  then  read  out. 

14.  The  input  and  fitted  spectra  are  plotted  in  a  graph  output. 

15.  The  above  process  is  repeated  for  each  run  in  the  data  set  until  the 
data  for  all  the  experimental  angles  has  been  processed. 

16.  When  the  preceding  analysis  is  completed  for  a  set  of  data  a  sub- 
routine entitled  LEASTSQ  may  be  called,  if  appropriate,  to  compute  the 
correlation  coefficients.   For  each  gamma  ray  in  the  spectrum  a  least 
squares  fit  of  the  counts  computed  at  the  experimental  angles  is  made  to 

the  usual  series-. 

1 
W(6)  =  2  A. P. (cos  0),   i  even 
i 
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where  P.  are  the  even  Legendre  polynomials.   The  fit  is  first  made  to 

An  +   A_P„  +  A. P.  then  A~  +  A„P„.   (Provision  has  been  made  for  fitting 
0    22    4  4       0    22 

to  the  series  terminating  with  P,  by  simply  changing  the  initial  value 

of  NP  to  4;  however,  see  reference  61  concerning  accuracy  of  P  terms). 

41 
The  computation,  which  follows  the  procedure  given  by  Rose   ,  proceeds 

as  follows: 

58 

a.   Table  gamma  ray  absorption  coefficients   in  copper  and  cor- 
responding energies  in  Mev.   Correct  intensities  and  sigmas 
for  attenuation  in  copper  target  backing. 

bo   Set  up  matrix  P. .  of  Legendre  polynomials  for  experimental 
angles. 

c.  Do  least  squares  calculation  to  find  coefficients  using  the 
reciprocal  of  the  estimated  intensity  variances  from  11  as 
weights. 

d.  Normalize  coefficients  and  correct  for  the  effect  of  the 
finite  solid  angle  subtended  by  the  scintillation  crystal. 

These  latter  corrections  are  precomputed  in  an  auxiliary 

41 
program  using  the  expressions  given  by  Rose.    Estimated 

variances  of  the  coefficient  ratios.,  A»/An,  etc.s  are  com- 

puted  from  the  expression  given  by  Reich  et  al.    (After 

2 

computation  variances  were  increased  by  a  factor  1/^   in 

2 
all  cases  where  the  computed  value  of  C  was  less  than  unity. ) 

e.  Calculate  anisotropy. 

2 

f.  Compute  the  criterion  £  .   This  quantitys  defined  by  Rose's 

analysis,  is  a  measure  of  non-statistical  errors  in  the  data. 

2 
If  such  errors  are  eliminated  the  value  of  C   serves  to  in- 
dicate whether  the  appropriate  number  of  polynomials  has  been 

used  in  the  correlation  function  expression.,  with  values  of 

2 
C   greater  than  unity  indicating  an  inappropriate  choice. 

It  is  pertinent  here  that  non-statistical  errors  have  already 
been  introduced  into  the  problem  in  11  above  so  that  this 
criterion  is  valid  only  to  the  extent  that  an  accurate  esti- 
mate of  the  true  variance  of  the  counts  has  been  achieved. 

2 
The  computed  values  of  £   in  general  are  near  unity  and  seem 

to  verify  that  the  estimated  variances  can  be  considered  to 

be  reasonably  good. 
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Selected  values  of  the  normalized  calculated  correlation 
function,  and  values  of  a  comparison  theoretical  function 
if  desired,  are  read  out,  and  similarly  normalized  experi- 
mental points  are  read  out.   (Experimental  values  are  un- 
corrected for  solid  angle,  the  correction  being  applied  to 
the  theoretical  function  for  comparison  with  experimental 
results  . ) 

Iterate  for  each  gamma-ray  energy. 

Reduce  number  of  polynomials  if  orders  higher  than  P  have 
been  used  and  iterate  solution,,   (The  program  also  calculates 
coefficients  for  the  cosine  series; 


m       _ 
W(9)  =  1/AQ  2  A2n  [COS  (0)] 


2n 


) 


jo   Return  to  main  routine . 

16 .  Repeat  entire  analysis  for  as  many  sets  of  data  as  indicated. 

17.  Combine  coefficients  to  get  weighted  mean  values  and  Chi  square 
values  if  more  than  one  set  of  data  is  analyzed.   Read  out  these  values 

18.  Combine  intensities  at  each  angle  from  all  sets  of  data  and  re- 
analyze for  coefficients  using  this  combined  data  set. 

C.   Glossary  of  names  of  selected  variables  and  constants  appearing 
in  program  SPECFIT.   An  asterisk  indicates  an  input  quantity. 


A(I,J) 

AINV(I,J) 
ANISO(I) 

ARGP(I) 
*B(I) 

*BACK(I) 

BACKGND(I) 

*BCHAN 


The  matric  product  of  the  response  matrix  with  its 
transpose,  weighted  with  the  WT  values. 

The  inverse  of  A(I,J). 

Anisotropy  for  Legendre  series  solution  (I  =  1) s 
and  for  cosine  series  solution  (I  -  2) . 

Vector  containing  cosines  of  the  angles  B3(I). 

A  vector  representing  the  product  of  the  inverse 
of  the  response  matrix  and  the  observed  spectrum. 

A  number  representing  background  counts  to  be 
subtracted  from  calibration  curve  I. 

The  calculated  background  spectrum  to  be  subtracted 
from  an  input  spectrum. 

The  lower  of  the  channel  pair  described  under  TCHAN 
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FIN(I,J,K) 
FMU(I) 

*HCM 
HOLD (I) 
LPKS 

*MAXCHAN 

*MESH 

*MINT 

*MINTB 

MPKS 

MPKS2 

*NBAOPT 


*NCUT 
NERFL 

*NHIEST 

*NLOWEST 

*NOCALB 
*NPHOTO 


Storage  for  FINTENS(I),  data  set  J3  run  K. 

Mass  absorption  coefficient  in  copper  for  EMU(I) 
Mev  gamma  ray  Units  in  cm^/gm. 

Distance  from  source  to  crystal  face  in  centimeters. 

A  temporary  storage  vector. 

Same  as  NPKS  but  not  including  any  straight  line 
component.  (See  NBAOPT) 

An  upper  channel  limit  in  loading  CALB  data.   Between 
MAXCHAN  and  MINT,  CALB(I)  values  are  filled  in  with 
the  value  of  CALB (MAXCHAN) . 

The  number  of  channels  in  input  spectrum. 

Number  of  channels  in  calibration  spectra. 

Number  of  background  channels. 

Initial  value  of  NPKS  for  restoring  the  value. 

Initial  value  of  LPKS  for  restoring  the  value. 

Background  option  selector. 

NBAOPT  -  0  No  background  read  in  nor  subtracted 
from  SPEC 

NBAOPT  =  1  Background  spectrum  is  read  in  and 

SPEC  is  corrected  for  this  background 

NBAOPT  =  2  Same  as  NBAOPT  =  1  and  in  addition  a 
straight,  line  is  included  in  the 
response  matrix. 

A  channel  number  below  which  the  graph  output 
values  are  divided  by  a  scale  factor  SCALE. 

An  output  of  matrix  inversion  subroutine.   A  value 

of  1  indicates  matrix  is  non  singular,  and  2  indicates 

a  singular  matrix. 

The  upper  channel  of  the  range  over  which  a  fit  is 
made  in  SPECFIT. 

The  lower  channel  of  the  range  over  which  a  fit  is 
made  in  SPECFIT. 

Number  of  calibration  spectra. 

The  photograph  number  of  the  observed  spectrum;  for 
identification  only. 
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*NPKS 


*NRUNS 
NP 

*NPRINT 


*NSETS 
OMEGA(IsJ) 

*PKSEP(I) 


*PHOT(I) 

*PKEN(I) 
POLY (I. J) 

R(I,J) 
RED 


*RES 

SAN (I) 
SAVE (I, J) 


The  number  of  spectral  components  with  which  the 
observed  spectrum  is  to  fitted.   Program  will  increase 
this  value  by  one  if  option  is  chosen  to  Include  a 
straight  line  in  R  (See  NBAOPT) . 

Number  of  runs  in  a  set  of  data. 

Number  of  Legendre  polynomials  to  be  fitted  in  first 
pass  through  LEASTSQ. 

Print  suppression  index 

NPRINT  =0  All  printouts  and  graphs  are  made 

NPRINT  =  1  Suppress  printouts  of  spectrum  values  and 
cosine  series  coefficients 

NPRINT  =  2  Same  as  1  and  also  suppress  graph  outputs 

NPRINT  =  3   Same  as  2  with  briefer  more  compact  arrange' 
ment  of  correlation  coefficients 

NPRINT  =  4  Same  as  3  but  graphs  plotted 
Number  of  sets  of  data. 

Used  for  temporary  storage  in  SPECFIT  and  diagonal 
elements  used  for  weights  in  LEASTSQ. 

The  number  of  channels  corresponding  to  1.022  Mev 
for  energy  calibration  of  CALB  curves.   Taken  as 
channel  distance  from  photo  peak  to  pair  peak  for 
higher  energy  calibration  gammas. 

The  channel  of  calibration  curve  I  in  which  the 
photopeak  occurs.   (May  be  an  interpolation  between 
two  integral  channels  numbers) . 

Energy  in  Mev  of  a  spectral  component. 

The  matrix  of  Legendre  polynomials  of  order  2(J-1) 
evaluated  at  angle  B3(I). 

The  response  matrix. 

A  normalizing  factor  for  3"  x  3"  unshielded  Nal(Tl) 
crystal  representing  the  probability  that  one  photon 
at  the  source  will  produce  one  count  in  the  photopeak 
channel. 

The  proton  energy  in  kev  at  the,  resonance  being 
investigated  (for  identification  only). 

Estimate  of  standard  deviation  of  ANISO(I) . 

An  auxiliary  matrix. 
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SBB(I) 

*SCALE 
SCOEFT(I) 
SIGCTS(I) 
SIGTENS(I) 

SIN(I,J,K,) 
*SPEC(I) 

STC(I3JSK) 

STORE (I, J) 

STSC(I9J,K) 

STW(ISJ) 

STSW(I,J) 

SUM7(I) 

*TBACK 

TC(I) 

*TCHAN 

TEMP (I) 

*TOPEN 

TS(I) 

*TX60SEC 

WT(I) 

XI(I) 


Coefficient  for  renormalizing  input  spectrum  with 
counts  from  an  auxiliary  monitor  (optional;  if  not 
entered  values  are  unity) . 

A  scale  factor  for  graph  output.   See  NCUT. 

The  standard  deviation  in  COEFT(I). 

The  calculated  error  in  COUNTS(I). 

The  error  for  calculated  spectral  intensity 
INTENS(I). 

Storage  for  errors  in  FIN(ISJ,K) 

The  vector  representing  the  counts  in  the  observed 
spectrum . 

Storage  for  accumulating  calculated  correlation 
coefficients  of  order  I,  for  gamma  ray  J  and 
data  set  K. 

A  temporary  storage  matrix  for  A  while  it  is  being 
inverted o 

Standard  deviation  of  STC(I,J,K)° 

Temporary  storage „ 

Temporary  storage „ 

Temporary  storage. 

Duration  of  background  run  in  units  of  seconds  x  60 

The  weighted  mean  value  of  correlation  coefficient 
of  order  I. 

The  upper  of  a  pair  of  input  spectrum  channels  to  be 

used  for  sensing  energy  calibration  of  the  input  spectrum, 

Temporary  storage. 

Energy  of  TCHAN  in  Mev. 

Standard  deviation  of  TC(I) 

Duration  of  a  run  in  units  of  seconds  x  60 „ 

The  weight  assigned  to  channel  I  in  least  squares 
solution  for  intensities. 

The  matric  product  of  POLY  and  Bl  weighted  with  the 
diagonal  elements  of  OMEGA. 
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XS(I)  A  vector  representing  the  energy  in  Mev  of  each 

channel  of  the  input  spectrum. 

Z(I)  A  temporary  storage  vector . 

*ZERO  The  number  which  if  any  diagonal  element  of  matrix 

A  is  smaller  than,  A  will  be  considered  singular . 
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D,   DATA  INPUT  SEQUENCE  FOR  PROGRAM  SPECFIT 


Input 

Number 
of  sets 

Number  of 
cards  in  set 

Input  name 

Remarks 

1 

1 

1 

NSETS,  NPRINT 

4 

2 

1 

1 

EL,  ECC,  HCM 

NOCALB,  MINT 


NOCALB 


MAXCHAN/7 


BACK(J)5  PHOT(J)3  ENCALB(J), 
MAXCHAN,  PKSEP(J) 


CALB(ISJ) 


NPKS,  NBAOPT 


6 

1 

NPKS 

PKEN(I) 

7 

1 

1 

BSLS  BENS  BPHOT,  MINTB, 
TBACK 

1 

8 

1 

MINTB/7 

BTGT(I) 

9 

1 

1 

NRUNSa  CTSNORA,  NLOWEST, 
NHIEST,  NCUTS  SCALE,  ZERO 

2 

10 

1 

1 

RES,  NPHOTO 

3 

11 

1 

1 

MESH,  TCHAN9  TOPEN,  BCHAN, 
BOTEN,  B3(NRUNDEX),  T60XSEC 

12 

1 

1 

CTSNORB 

13 

1 

MESH/ 7 

SPEC(I) 

14 

1 

4 

Graph  data  cards 

5 

REMARKS: 

1.)   Read  in  only  for  NBAOPT  greater  than  1„ 

2o)   Input  resumes  here  for  each  successive  data  set. 

3.)   Input  resumes  here  for  each  successive  run  within  a  data  set. 

4.)   All  loading  formats  are  multiples  of  F10»3  and  110  except  that  ZERO 
is  ElOol  and  TBACK  is  F20.3o 

5o)   Blank  cards  may  be  loaded  if  graph  option  not  used. 
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E.   Description  of  program  CONTOUR 

As  described  in  appendix  1  the  correlation  function  coefficients  for 
(p,<f  )  and  (p,xtf)  intermediate  unobserved  reactions  for  spin  1/2  target 
nuclides  can  be  described  as  a  product  of,  at  most,  three  functions  of  the 
form: 


+  bik<Sk +  cik£ 


8ik  -  "  ,       rl 

1  +dk 


where   i  denotes  the  order  of  the  coefficient  and  k  is  the  index  of 
the  mixing  parameter  £   .   For  computation  the  substitution  O  =  tan  Y 
is  commonly  used,  and  leads   to  g  =  a  +  (3  sin  (2Y  +T) 


.  ,  a  +  c 

with     a   =  — -x — 


P  =  +  |  [b2  +  (a  -  c)2] 

JT  =  tan"1  (±r=-£  ) 
b/  |b| 


1/2 


s  = 


where  indices  have  been  dropped  for  clarity. 

This  program  computes  sets  of  possible  values  of  two  mixing  para- 
meters for  various  fixed  values  of  correlation  coefficients.   Using  a 
graph  plotter  output  medium  available  through  a  locally  developed  sub- 
routine, these  sets  of  values  are  plotted  to  form  a  set  of  contours 
along  which  the  value  of  the  correlation  coefficient  is  constant. 
Steps  in  the  computation  are  as  follows: 

1.  Read  in  data  for  one  possible  spin  sequence  to  be  examined . 
Compute  values  of  (X,    p  and  #   for  the  lowest  order  coefficient 
which  occurs.   If  three  mixing  parameters  occur  simultaneously, 
read  out  the  expression  for  the  fixed  g  as  a  function  of  the 
fixed  mixing  parameter.   (All  loading  formats  are  standardized 
to  be  F10.3  and  110  for  floating  and  fixed  point  variables 
respectively. ) 

2.  Find  and  read  out  the  maximum  and  minimum  possible  values  of 
the  coefficient  and  the  corresponding  values  of  the  Y's. 
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Fix  the  coefficient  at  an  initial  value  within  its  range. 
Assign  values  to  one  of  the  Y's  in  equal  increments  through- 
out the  range  =90   to  +90  and  solve  for  the  corresponding 
values  of  the  remaining  Y,  which  in  general  will  be  double 
valued  over  the  region  selected. 

Read  out  pairs  of  Y  values  for  the  selected  coefficient 
value  both  to  the  print  medium  and  to  the  graph  plotter. 
Repeat  steps  three  and  four  for  successively  larger  values 
of  the  coefficient  until  a  specified  number  of  contours  has 
been  plotted. 

Repeat  steps  two  through  five  until  contour  plots  for  co- 
efficients of  all  orders  permitted  by  the  given  spin  sequence 
have  been  produced . 


F.   Input  variables  for  CONTOUR 
NSETS 


Number  of  spin  sequences  for  which  plots  are  to  be 
made . 


JFIND 


JSTEP 
JFIX 

IMAX 

JMAX 

NCONTS 
NPTS 


Index  of  the  mixing  parameter  to  be  solved  for.  (Y 
values  on  graph) . 

1  =  proton  partial  wave  mixing 

2  =  unobserved  gamma  multipolarity  mixing 

3  =  observed  gamma  multipolarity  mixing 

4  =  channel  spin  mixing 

Index  of  mixing  parameter  assigned  successive  values 
(x  values  on  graph)  See  JFIND „ 

Index  of  mixing  parameter  to  be  held  constant  when 
three  parameters  occur  simultaneously.   (Set  to  five 
for  (p, $  )  reaction  where  only  two  parameters  occur.) 

Index  specifying  the  highest  order  coefficient  occurring 
in  the  spin  sequence.  (Order  =  2»IMAX) 

Largest  value  which  the  index  of  mixing  ratio  assumes 
for  use  in  reading  input  data.  (Usually  set  at  four.) 

Number  of  contours  to  be  plotted  on  each  graph. 

Number  of  pairs  of  points  to  be  computed  for  each  branch 
of  one  contour. 


DELTA (JFIX) 


Fixed  value  assigned  to  the  third  mixing  parameter  during 
a  calculation. 
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RHO  Energy  dependent  Coulomb  phase  factor  to  be  multiplied 

into  the  b  coefficients  for  proton-partial-wave 
mixing. 

A(I,J)  The  a,  b,  and  c  factors  occurring  in  the  expression  for 

B(I,J)  go   I  denotes  the  coefficient  order  (order  -  21)  and 

C(I9J)  J  the  mixing  parameter  index.   See  FIND,, 
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